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An  LKB  10700-2  batch  microcalorimeter  was  modified  to  improve  its  accuracy  and 
reproducibility  and  to  shorten  the  time  required  for  each  experiment.  The  excess  enthalpies, 
partial  molal  volumes  and  partial  molal  expansibilities  of  aqueous  solutions  of  N-methylformamide, 
N-methylacetamide,  N-methylpropionamide,  N-butylacetamide,  ethylene  glycol,  pentaerythritol, 
glucose,  sucrose,  urea  and  binary  mixtures  of  these  compounds  have  been  measured  at  25°C.  A 
statistical  analysis  of  the  pairwise  heats  of  interaction  calculated  from  the  excess  enthalpies  showed 
that  these  interaction  heats  could  be  roughly  represented  by  a sum  of  several  independent  group 
interactions.  The  sign  and  magnitude  of  the  following  group  interactions  are  reported: 


O 


O 


(CH2)-(CH2),  (C-NH)-(C-NH),  (CHOH)-(CHOH), 

o O 

II  II 

(CH2)-(CHOH),  (CH2)-(C-NH)  and  (C-NH)-(CHOH). 


These  group  interaction  values  were  used  to  calculate  the  pairwise  heats  of  interaction  for  several 
previously  reported  nonelectrolytes  with  a a of  45  cal  kg/mol^. 
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THE  ENTHALPY  OF  INTERACTION  OF  NONELECTROLYTE 
PAIRS  IN  DILUTE  AQUEOUS  SOLUTIONS 

I.  INTRODUCTION. 

Many  theories  exist  to  explain  the  behavior  of  ions  in  aqueous  as  well  as  nonaqueous 
solutions;  yet  the  development  of  such  theories  to  explain  the  behavior  of  nonelectrolyte  solutions 
is  just  beginning.  This  is  in  spite  of  the  fact  that  the  behavior  of  nonelectrolyte  solutions  should  be 
easier  to  interpret  since  long-range  forces  due  to  charges  are  absent. 

The  overall  goal  of  this  research  was  to  understand  how  individual  functional  groups 
affect  the  interactions  between  pairs  of  nonelectrolytes  in  aqueous  solutions.  Around  each 
nonelectroly  te  molecule  there  will  be  one  or  more-layers  of  water  molecules  whichdifferfrom  the 
bulk  water;  this  is  the  cosphere  model.1  When  the  cosphere  on  one  nonelectrolyte  molecule 
overlaps  the  cosphere  on  another  molecule,  the  resulting  interaction  forces,  attractive  or  repulsive, 
were  expected  to  depend  on  dipole-dipole  effects,  van  der  Waal  forces,  the  effect  due  to  the 
replacement  of  water  hydrogen  bonds  with  intermolecular  hydrogen  bonds  and  cosphere  overlap 
forces.1  In  this  work  the  heats  of  interaction  resulting  from  a dilution  of  a nonelectrolyte  solution 
were  measured.  The  heats  were  measured  since  very  little  is  known  about  the  interactions  between 
nonelectrolytes  in  aqueous  solutions  and  because  enthalpy  measurements  are  easily  performed  with 
the  equipment  available. 

Kauzmann  demonstrated2  the  tendency  of  nonpolar  groups  in  proteins  to  adhere  to  one 
another  in  aqueous  environments.  This  is  termed  hydrophobic  bonding.  Qualitative  results  found 
in  the  work  of  Scheraga3  and  the  free  energy  data  for  the  large  tetraalkylammonium  salt  systems4 
suggests  that  hydrophobic  bonds  are  attractive.  Based  on  the  hydrogen  bonding  found  in  amide 
solutions5  '9  one  would  expect  an  attractive  interaction  between  these  groups.  In  this  study  we  are 
trying  to  measure  quantitative  values  for  the  heats  of  these  interactions.  Quantitative  values  for  the 
heats  of  interactions  between  different  molecular  groups  should  provide  accurate  predictions  for 
the  thermodynamic  behavior  of  large  molecules  that  are  composed  of  pieces  for  which  the 
interactions  are  known.  A large  variety  of  molecules  could  conceivably  be  predicted  with  a small 
amount  of  data.  For  example,  once  the  interaction  between  a hydrocarbon  with  a hydrocarbon,  a 
hydrocarbon  with  a (CHOH)  group  and  a (CHOH)  with  a (CHOH)  group  are  known  it  should  be 
possible  to  predict  the  properties  of  a wide  variety  of  alcohols,  hydrocarbons  and  carbohydrate 
mixtures  simply  on  the  basis  of  these  three  parameters.  When  the  interaction  parameters 
for  a wide  variety  of  functional  groups  are  measured,  it  should  be  possible  to  understand  how 
various  nonelectrolytes  interact  with  proteins  to  alter  their  relative  stability  in  aqueous  solutions. 

It  may  even  be  possible  to  predict  for  example  the  ability  of  an  enzyme  to  bind  to  a substrate  and 
perhaps  predict  small  changes  in  the  known  enzyme  structures  for  this  process.  This  research  will 
not  achieve  such  far-ranging  goals,  but  these  objectives  are  possible  and  as  a start,  the  present 
research  has  collected  data  on  a wide  variety  of  nonelectrolytes  to  determine  if  accurate  predictions 
of  the  heats  of  interaction  on  nonelectrolytes  can  be  made  from  a knowledge  of  the  individual 
group  interactions  involved  by  using  a simple  additivity  rule. 

Thermodynamic  investigations  of  nonelectrolyte  systems  have  been  reported  by 
several  authors,  Kozak.  Knight  and  Kauzmann,10  Friedman  and  Krishnan.1  1 Lilley  and  Scott,* 
Lange  et  a/.,12’13  Franks,**  Wood  et  a/.,1415  and  others.  Much  of  the  earlier  data  are  not 

* Lilley,  T.  H.,  and  Scott,  R.  H.  Private  communication,  1974. 

**  Franks,  F.  Private  communication,  1975. 
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useful  in  investigating  the  heats  of  interaction  between  nonelectrolyte  pairs  and  their  mixtures  in 
dilute  solutions.  Some  of  the  data  have  not  been  measured  with  enough  accuracy  or  have  been 
obtained  with  equipment  not  suitable  for  heats  of  dilution.  Most  of  the  enthalpy  data  consists  of 
dilutions  of  single  nonelectrolytes.  The  aim  of  this  research  was  to  measure  the  effect  that  different 
functional  groups  have  on  the  observed  pairwise  heats  of  interaction.  This  could  only  be 
accomplished  by  testing  calorimetrically  dilutions  of  binary  nonelectrolyte  mixtures  for  which  the 
data  is  extremely  scarce.  Also  it  is  only  by  a large  systematic  study  such  as  performed  here  that 
the  formulation  of  overall  trends  in  the  data  and  quantitative  results  could  be  expected. 

The  compounds  chosen  for  this  study  were:  N-methylformamide  (NMF),  N-methyl- 
acetamide  (NMA),  N-methylpropionamide  (NMP),  N-butylacetamide  (NBA),  urea  (U),  ethylene 
glycol  (EG),  pentaerythritol  (P),  glucose  (G),  and  sucrose  (S).  Aqueous  solutions  of  these 
compounds  as  well  as  tlfeirT-Pmlxfutes  were  diluted  from  2 molal  to  approximately  0.2  molal  or 
lower.  McMillan-Mayer  theory  1 6 was  used  to  show  that  measurements  performed  in  this  way 
resulted  in  the  heats  of  interactions  for  pairs  of  solute  molecules.  These  heats  of  interaction  were 
then  analyzed  by  a general  least  squares  program  to  determine  the  quantitative  values  for  the 
pairwise  heats  of  interaction  between  hydrocarbon  groups,  amide  groups,  (CHOH)  groups  and  the 
unlike  group  interactions.  With  these  data  we  were  able  to  obtain  for  the  first  time  a quantitative 
prediction  of  the  enthalpies  of  interactions  for  a number  of  simple  nonelectrolytes  containing 
different  amounts  of  hydrocarbons,  amide  and  (CHOH)  groups. 

11.  EXPERIMENTATION. 


The  experimental  procedure  and  results  of  the  density  measurements  are  presented  in 
this  section  as  well  as  a description  of  the  equipment  used  to  measure  the  heats  of  dilution.  The 
procedure  and  method  of  calculation  of  the  heats  of  dilution  are  also  presented. 

A.  Solutions. 

Nine  compounds  were  tested  for  the  heat  of  dilution  and  density  experiments.  These 
were  the  following:  four  amides  - N-methylformamide  (NMF),  N-methylacetamide  (NMA). 
N-methylpropionamide  (NMP),  and  N-butylacetamide  (NBA);  four  polyhydroxy  compounds  - 
ethylene  glycol  (EG),  pentaerythritol  (P),  glucose  (G),  and  sucrose  (S);  and  urea  (U).  The  NMP  and 
NBA  were  Eastman  Kodak  high-purity  compounds  and  the  NMF  was  “Baker”  grade  material;  these 
compounds  were  used  as  received.  The  NMA  was  Aldrich  analytical  grade  which  was  recrystallized 
three  times  in  a dry  box  and  then  zone-melted  for  five  days  before  use.  Ethylene  glycol  was 
Fischer-certified  and  used  as  received.  Glucose,  sucrose  and  urea  were  Fischer-certified  and  dried 
in  a vacuum  oven  for  three  days  at  60°C  before  use.  The  pentaerythritol  was  Baker  analytical  grade 
which  was  recrystallized  from  water  three  times  and  dried  in  a vacuum  oven  for  three  days  at  90°C 
before  use.  Table  1 lists  the  water  content  of  all  the  compounds,  as  determined  by  Karl  Fischer 
analysis,  and  table  2 the  impurities  of  the  Fischer  certified  compounds  as  listed  on  the  bottle  labels. 
The  NMF,  NMP  and  NBA  were  analyzed  (gas-chromatography)  and  found  to  be  greater  than  99.3% 
pure.  The  only  impurity  found  was  water  which  was  corrected  for  in  making  up  these  solutions. 
The  pentaerythritol  which  had  a melting  range  of  252-258°C  before  purification  melted  at  260°C 
after  purification. 


Table  1 . Water  Content  by  Karl  Fischer  Analysis 


Compound 

Wt  % h2o 

N-Methylformamide 

0.50 

N-methylacetamide 

<0.03 

N-Methylprapionamide 

0.60 

N-Butylacetamide 

0.19 

Urea 

0.03 

Ethylene  glycol 

0.06 

Pentaerythritol 

0.06 

Glucose 

0.08 

Sucrose 

<0.03 

Table  2.  Fischer-Certified  Label  Impurities 


Compound 

Urea 

Ethylene  glycol 

Glucose 

Sucrose 

Insoluble  matter 

0.009* 

0.002 

0.003 

Acidity 

0.009 

0.01 

0.002 

Chloride 

0.003 

0.003 

0.003 

Sulfate  and  sulfite 

0.001 

0.002 

0.001 

Arsenic 

0.00001 

Heavy  metals  (Pb) 

0.0003 

0.0002 

0.002 

Iron 

0.0002 

0.0001 

0.0001 

Invert  sugar 

0.05 

• % Impurity. 
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All  solutions  of  these  compounds  were  prepared  by  weight  to  ±0.01%  with  distilled 
water  which  was  trickled  through  a Barnstead  ion  exchange  column.  Bouyancy  corrections  were 
made  in  all  work.  In  all  cases  the  solutions  were  made  up  and  used  within  a 36-hr  period  during 
which  time  they  were  stored  in  polyethylene  screw  cap  bottles.  This  treatment  of  the  solutions 
was  thought  advisable  to  forestall  any  decomposition  which  might  occur  due  to  bacterial 
contamination.  This  is  especially  important  in  the  case  ot  the  sugars.  In  three  experiments,  a 
solution  of  sucrose  was  measured  approximately  52  hours  after  preparation.  The  heat  of  dilution 
was  an  average  of  2%  low. 

B.  Density  Measurements. 

The  solution  densities  of  1 and  0.3  moial,  of  the  compounds  used  in  this  study,  were 
measured  at  0,  25  and  45°C  using  a single-stem  dilatometer.  The  density  of  pentaerythritol  (P)  was 
measured  at  0,  25  and  45°C  for  only  one  concentration  (0.3  ni)  because  of  its  low  solubility  in 
water.  The  dilatometer  was  constructed  from  a 25-cc  volumetric  flask  to  which  a 21-cm-long  stem 
of  precision  bore  tubing  (approximately  1-mm  inside  diameter)  was  attached.  The  dilatometer  was 
loaded  with  degassed  material  by  a syringe  equipped  with  Teflon  tubing.  The  experiment  measures 
volumes  of  a known  weight  of  solution  using  the  following  general  equation: 

V = V°  - A X L (1) 


where 


V = volume  of  material  in  cc  of  a known  weight  in  gm 
V°  = volume  of  dilatometer  to  top  of  stem  in  cc 
A = stem  area  in  cm  not  filled  with  material 
L = stem  length  in  cm  not  filled  with  material. 

Calibrations  were  performed  at  0.  25  and  45°C  on  degassed,  ion-exchanged,  distilled  water  by  use 
of  published  density  data.1 7 Bath  temperatures  were  monitored  within  ±0.02°C  with  a platinum 
resistance  thermometer.  The  calibrations  provided  the  capillary  reading  L as  a linear  function  of  the 
total  volume  shown  in  the  following  equations  for  the  three  temperatures: 

0°C  V = 25.3956(3)  - 0.01833  X L (2) 

25°C  V = 25.4016(5)  - 0.01828  X L (3) 

45°C  V = 25.4079(5)  - 0.01826  X L (4) 

w'here  V and  L have  been  described  in  equation  1 and  the  numbers  in  the  parentheses  are  the 

standard  deviations  of  the  last  digit  for  a single  measurement,  i.e.,  25.3956(3)  = 25.3959  to 
25.3953. 


The  solutions  were  prepared  by  weight  using  deaerated  water.  The  dilatometer  was 
rinsed  several  times  with  water  followed  by  Fischer-certified  acetone  and  dried  inside  with  a stream 
ot  dry  nitrogen.  The  outside  of  the  dilatometer  was  then  wiped  with  a slightly  damp  towel  and 
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placed  in  the  balance  case  to  equilibrate  for  approximately  one  hour.  The  empty  weight  of  the 
dilatometer  was  recorded  when  two  successive  weighings,  five  minutes  apart,  agreed.  The  apparatus 
was  then  filled  with  solution,  again  allowed  to  equilibrated  in  the  balance  case  and  the  full  weight  of 
the  dilatometer  recorded  when  two  successive  weighings  agreed.  Measurements  of  stem  length  not 
filled  with  solution  (L)  were  then  made  at  0,  25  and  45°  ±0.01°C  with  a caliper  gauge  to  ±0.01  cm. 

Thermal  equilibrium  of  the  dilatometer  and  its  contents  with  each  bath  temperature  was  achieved 
when  three  successive  height  measurements  performed  five  minutes  apart  agreed.  Because  the 
expansibility  of  the  solutions  exceeded  the  capacity  of  the  dilatometer  over  the  measured 
temperature  range,  all  the  45°C  height  measurements  were  made  after  some  solution  had  overflowed 
from  the  dilatometer.  The  outside  of  the  dilatometer  was  then  rinsed  with  distilled  water,  wiped 
with  a damp  towel,  allowed  to  equilibrate  in  the  balance  case  and  weighed.  The  solution  densities 
were  calculated  by  dividing  the  solution  weight  by  the  measured  solution  volume.  The  estimated 
precision  of  the  density  measured  was  1 X 10*^  gm/cm^. 

The  apparent  molal  volume  of  the  1 and  0.3  molal  solutions  were  calculated  at  0,  25  and 
45°C  from  the  solution  densities  using  the  following  equation: 
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where 


0V  = apparent  molal  volume  in  cm^/mol 

d = solution  density  at  desired  temperature  in  gm/cm^ 

MW  = molecular  weight  of  solute  in  gm 

m = molality  of  solution  in  mol/kg 

d°  = density  of  water  at  desired  temperature  in  gm/cm^. 

An  extrapolation  to  infinite  dilution  of  0V  versus  m plots  at  0,  25  and  45°C  was  used  to  obtain  the 
partial  molal  volumes  at  infinite  dilution  <t>y  for  the  compounds  measured.  These  values  of  <p°  were 
then  fitted  to  a second  order  polynomial  of  the  form: 

</>v  = a + bt  + ct^ 

where  a,  b and  c are  constants  and^t  is  the  temperature  in  degrees  centigrade.  The  partial  molal 
expansibility  at  infinite  dilution  </>g,  was  then  calculated  for  each  solute  by  differentiating  these 
polynomials  with  respect  to  temperature.  The  results  as  well  as  available  literature  data  are 
presented  in  table  3. 

Of  the  nine  compounds  studied,  N-methylformamide  (NMF)  and  N-butylacetamide  (NBA) 
were  the  only  compounds  for  which  literature  density  data  could  not  be  found.  The  data  are  shown 
in  table  3 are  ±0.02  based  on  the  precision  of  the  density  measurement.  Comparison  of  the  present 
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f Table  3.  The  Partial  Molal  Volume  (0V)  and  Partial  Molal  Expansibility  (0g) 

of  Aqueous  Nonelectrolyte  Solutions 


Compound 

m 

■ 

0£ 

■ 

mol/kg 

cm  /mol 

o°c 

25°C 

45°C 

0°C 

25°C 

45°C 

N-Methylformamide  (NMF) 

54.8 

56.6 

57.8 

54.9 

56.6 

57.9 

0 expt. 

54.9  ( 3)* 

56.6  ( 3) 

57.9  ( 3) 

.070 

.066 

.064 

N-Methylacetamide  (NMA) 

1.0003 

73.3 

74.5 

1.0000 

71.9 

0.3000 

72.5 

73.5 

74.7 

0 expt. 

72.7  ( 3) 

73.6  ( 3) 

74.8  ( 3) 

.023 

.049 

.070 

0 lit.1  8 

73.9  ( 1) 

0 lit.1  9 

72.9  ( 4) 

N-Methylpropionamide  (NMP) 

0.9977 

86.8 

88.7 

90.3 

0.2981 

86.0 

89.3 

90.8 

0 expt. 

88.4  ( 3) 

89.5  ( 3) 

90.9 

.030 

.058 

.082 

0 lit  7° 

90.4  ( l) 

N-Butylacetamide  (NBA) 

1.0001 

117.4 

120.9 

123.9 

0.2987 

119.1 

121.6 

124.2 

0 expt. 

119.8  ( 3) 

121.9  ( 3) 

124.3  ( 3) 

.064 

.104 

.136 

Urea  (U) 

1 .0000 

42.1 

44.3 

45.4 

0.3001 

41.6 

43.9 

45.1 

0 expt 

41.4  ( 3) 

43.8  ( 3) 

44.9  ( 3) 

.119 

.073 

.037 

0 lit.21 

41.75  (10) 

44.20(10) 

45.10(10) 

.116 

.073 

.039 

0 lit.22 

41.9  ( 3) 

44.2  ( 3) 

45.3  ( 3) 

.113 

.072 

.039 

0 lit.1  8 

43.8  ( 1) 

0 lit.23-24 

43.7  ( 3) 

Ethylene  glycol  (EG) 

1.0111 

53.3 

54.5 

55.4 

0.9993 

53.4 

54.5 

0.3003 

53.4 

54.6 

55.5 

0 expt. 

53.4  ( 3) 

54.6  ( 3) 

55.6  ( 3) 

.047 

.049 

.051 

0 lit. 25 

54.2  ( 2) 

0 lit.26 

54.3  ( 3) 

Pentaerythritol  (?) 

0.3000 

99.7  ( 3) 

102.0  ( 3) 

103.6  ( 3) 

.099 

.085 

.075 

0 lit.27 

101.3  ( 5) 

0 lit.28 

101.9  ( 1) 

0 lit.29 

101.9  ( 1) 

Glucose  (G) 

0.9869 

109.4 

112.3 

113.9 

0.3005 

109.0 

112.2 

114.0 

0 expt. 

108.8  ( 3) 

112.2  ( 3) 

114.0  ( 3) 

.162 

.110 

.070 

0 lit.  2 

108.8  ( 1) 

112.1  ( 1) 

113.7  ( 1) 

.161 

.103 

.057 

Sucrose  (S) 

1.0002 

208.1 

212.6 

215.2 

0.3000 

206.7 

211.5 

214.5 

0 expt. 

206.3  ( 3) 

211.0  ( 3) 

214.2  ( 3) 

.204 

.173 

.148 

0 lit*0 

206.4  ( 1) 

211.5  ( 1) 

0 lit.22 

206.3  ( 1) 

211.0  ( 1) 

213.6  ( 1) 

.220 

.156 

.104 

* The  number  in  the  parenthesis  represents  the  estimated  error  of  the  last  digit,  eg.,  54.9  (3)  = 54.6  to  55.2. 
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data  tor  urea,  sucrose,  and  glucose  with  literature  data  at  0°C  and  25°C  demonstrates  that  this 
accuracy  was  achieved.  The  data  at  45°C  appear  to  be  considerably  less  accurate.  This  is 
probably  a result  of  the  experimental  procedure  used  to  measure  the  density  at  45°C.  The 
expansibility  of  all  the  solutions  exceeded  the  capacity  of  the  dilatometer  at  this  temperature  and 
its  contents  after  a measurement  at  45°C.  The  errors  involved  in  this  procedure  would  account  for 
the  reduced  accuracy  noted  in  for  the  45°data. 

Hepler3 1 and  Neal  and  Goring22  have  reported  that  the  <p £ of  a substance  in  water  could 
be  used  as  an  index  of  its  structural  interaction  with  water.  They  considered  water  to  be  an 
equilibrium  mixture  of  bulky  icelike  species  with  a relatively  low  density  and  a denser  species  formed 
by  breaking  or  deforming  the  hydrogen  bond  network  of  the  bulky  species.  Hydrophobic  solutes 
are  believed  to  cause  an  increase  in  water  structure  around  the  solute  molecules;  when  the  solution 
temperature  is  increased,  the  change  in  <t> ^ would  be  less  than  expected  for  an  ideal  solute.  The 
opposite  behavior  is  expected  for  hydrophilic  solutes  which  break  water  structure  around  the 
solute.  This  produces  larger  than  expected  values  in  0£. 

Neal  and  Goring  believe  that  the  magnitude  of  (the  apparent  specific  expansibility), 
shows  good  correlations  with  the  ability  of  a solute  to  alter  water  structure.  Large  values  of 
were  used  to  identify  solutes  as  water  structure  breakers  while  smaller  or  negative  values  indicated 
structure-making  solutes. 

Hepler3 1 used  the  results  of  several  literature  studies  on  the  density  of  aqueous  solutions 
to  show  that  the  sign  of  (6</>£/5T)p  is  positive  for  structure-making  solutes  and  negative  for 
structure  breakers.  He  also  considered  the  magnitude  of  this  slope  term  to  be  indicative  of  the 
ability  of  a solute  to  make  or  break  water  structure. 

Based  on  these  theories,  glucose,  sucrose  and  urea  were  classified  as  structure  breakers 
because  each  of  these  solutes  had  large  values  for  0£o  and  exhibited  a fairly  large  negative 

O O 

change  in  0£  with  respect  to  temperature.  The  values  of  (50g/6T)p  for  N-methylacetamide. 
N-methylpropionamide  and  N-butylacetamide,  shown  in  table  3,  are  fairly  large  and  positive.  These 
solutes  can  therefore  be  classified  as  structure  makers.  The  results  of  other  workers  using  nuclear 
magnetic  resonance  spectroscopy3  2 and  dielectric  constant  measurements3  3 have  also  shown  that 
these  amides  appear  to  make  water  structure.  N-Methylformamide  and  ethylene  glycol  do  not 
show  a sufficient  0£  slope  change  to  classify  these  solutes  as  structure  makers  or  breakers  whereas 
pentaerythritol  may  be  a structure  breaker  due  to  its  positive  0£  slope.  The  magnitude  of  this 
slope  change  however  is  much  smaller  than  that  found  for  other  structure-breaking  solutes  such  as 
urea  and  glucose. 

The  structure  theories  of  Neal  and  Goring  and  Hepler  encounter  problems  when  the 
results  obtained  for  ethylene  glycol  and  pentaerythritol  are  examined  more  closely.  Hydrophilic  OH 
groups  in  a molecule  were  considered  to  have  a structure-breaking  effect  by  these  authors. 
Pentaerythritol  and  ethylene  glycol  should  therefore  exhibit  structure-breaking  behavior  because 
of  the  number  of  hydrophilic  OH  groups  in  these  molecules.  Yet  based  on  Neal’s  and  Goring’s  and 
Hepler’s  theory,  ethylene  glycol  is  neither  a structure  breaker  or  maker  and  pentaerythritol  is  only 
a marginal  structure  breaker.  The  validity  of  classifying  amides  as  structure  makers  is  also 
questionable  since  other  experimental  techniques  have  been  used  to  classify  these  compounds  as 
either  water  structure  makers3  2 or  breakers.34 
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Finally,  even  the  validity  of  using  data  to  determine  a solutes  effect  on  water  structure 
is  uncertain  since  Holtzer  and  Emerson34  have  questioned  the  usefulness  of  water  structure  itself  in 
rationalizing  the  properties  of  aqueous  solutions. 

C.  The  Calorimeter. 

All  heat  of  dilution  experiments  were  performed  using  an  LKB  10700-2  batch 
microcalorimeter  previously  described  by  Wadso,36  LKB  Instruments’7  and  Levin.38  It  is  of  the 
heat  flow  design  utilizing  a rotatable  calorimeter  unit  suspended  in  a thermostated  air  bath.  The 
rotatable  unit  consists  of  an  aluminum  block  serving  as  a heat  sink.  Block  temperature  is 
continuously  monitored  with  a thermistor  bridge.  Two  18-carat-gold  reaction  cells  are  located  in 
the  aluminum  block  with  each  cell  divided  into  two  compartments  whose  volumes  are  approximately 
3 and  5 cc  respectively.  During  an  experiment,  the  rotation  of  the  calorimeter  unit  allows  the 
contents  of  the  cell  compartments  to  mix.  Heaters  located  within  each  cell  are  used  to  calibrate 
the  calorimeter.  Each  reaction  cell  is  sandwiched  between  two  semiconductor  thermopiles.  The 
difference  in  temperature  produced  between  the  sample  and  reference  cells  during  an  experiment 
causes  a voltage  difference  which  is  proportional  to  the  difference  in  heat  flow  between  the  cells 
through  the  thermopiles  to  the  aluminum  block.  The  heat  of  reaction  is  proportional  to  the 
integral  of  the  heat  flow.  The  thermopiles  sense  this  temperature  difference  as  a voltage  which  is 
then  amplified  by  a Keithley  instruments  model  150B  microvoltmeter  and  fed  into  a Sargent 
recorder,  equipped  with  a ball-and-disk  integrator.  An  example  of  a typical  heat  of  dilution  curve 
is  given  in  figure  1 . The  area  under  this  curve,  when  multiplied  by  the  calorimeter  constant  is  a 
direct  measurement  of  the  heat  produced  in  an  experiment. 

Since  the  calorimeter  was  not  used  for  approximately  one  year,  it  was  thoroughly  checked 
for  proper  operation  before  the  start  of  any  experiments.  The  thermostated  air  bath,  which  was 
inoperative  due  to  a broken  heater, and  cooler,  was  modified  by  using  a drinking  fountain  to  supply 
the  calorimeter  with  1 5°C  cooling  water.  This  change  was  made  because  no  cost  was  involved  for 
the  cooler  and  any  repairs  needed  could  be  easily  made  by  the  maintenance  staff  of  the  University. 
The  85-ohm  nichrome  heater  for  the  air  bath  was  found  to  be  shorted  out  and  was  modified  in 
design  to  prevent  possible  future  electrical  shorts.  A voltmeter  was  inserted  into  the  air  bath  heater 
circuit  to  allow  continuous  monitoring  of  the  amount  of  power  supplied  to  the  air  bath  heater. 
Previously  a panel  light  was  used  to  indicate  whether  the  heater  was  full  on  or  off.  An  adjustable 
resistor  (Rmax  = 100  ohms,  Rexp  = 10  ohms)  was  mounted  externally  and  wired  in  series  with  the 
85-ohm  air-bath  heater.  As  a result  of  this  change  it  is  now  possible  to  adjust  the  85-ohm  air-bath 
heater  to  half  power  using  the  external  resistor  and  continuously  monitor  the  power  level  with  the 
added  voltmeter.  Coolant  water  temperature  was  set  to  provide  an  exit  water  temperature  of  1 5°C 
with  a flow  rate  of  about  0.2  liters  per  minute  for  25°C  bath  operation  (860  on  thermistor  dial). 

The  air  bath  and  block  temperatures  were  measured  at  the  above  setting  with  a thermocouple, 
using  an  ice  bath  reference  junction,  checked  against  a platinum  resistance  thermometer  in  a 
25.00-  ± -0.01°-C  water  bath.  The  thermocouple  voltage  read  with  the  Keithley  microvoltmeter 
indicated  the  calorimeter’s  air  bath  and  block  temperature  was  25.00  ±0.02°C.  Another  check  of 
this  temperature  was  provided  with  a microvoltmeter  null  reading  when  the  thermocouple  reference 
junction  was  placed  in  the  air  bath  while  the  primary  junction  was  mounted  in  the  block.  The 
resistance  of  the  calibration  heater  located  in  the  sample  and  reference  cells  was  measured  with  an 
ESI  voltohmmeter  and  found  to  be  50.060  ohms  and  49.827  ohms  respectively.  These  measured 
values  compared  favorably  with  the  previously  measured  values  of  50.06  ohms  and  49.835  ohms  by 
Dr.  R.  Wood  on  20  September  1971 . The  calorimeter  has  a power  supply  which  provides  constant 
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current  inputs  to  the  cell  calibration  heaters.  In  the  past  these  current  values  and  resistances  of  the 
calibration  heaters  would  be  used  with  the  areas  of  the  calibration  curves  to  calculate  calorimeter 
calibration  constants  for  all  experiments  from  the  following  equation: 

(I)2  (R)  (0 

- — ■ ■■  (o) 

eR  (4. 184)  (A) 

where 

eR  = calibration  constant  for  a particular  Keithley  range  given  in 
calories  per  Keithley  output  (Vs) 

I = current  values  in  ma  of  the  current  source.  Values  measured 
before  the  start  of  experiments  and  assumed  constant 

R = resistance  in  ohms  of  calibration  heater  in  sample  cell 

t = time  in  seconds  current  applied  to  cell  heater 

A = area  in  volt  seconds  under  heat  curve  produced  by  electrical 
heating 

These  calibration  constants  would  then  be  multiplied  by  the  area  under  a dilution  curve  to  give  the 
heat  of  dilution  in  calories.  There  existed  no  easy  method  of  checking  these  current  values  on  a 
daily  basis.  Binding  posts  were  connected  into  the  constant  current  circuit  to  allow  voltage 
monitoring  with  an  ESI  voltohmmeter  whenever  an  electrical  current  was  input  to  a cell  heater 
during  a calibration.  Because  of  the  small  lead  resistance  between  the  binding-post  connection 
and  the  calorimeter  heater,  the  following  equation  was  used  for  calculating  heats: 


(VESl)2 

CR ; (7) 

(R+Rleads)2  (A)  (4. 184) 

where 

R,  t and  A have  been  described  in  equation  6 
^ESI  = voRage  at  binding  posts  measured  with  ESI 

R+Rleads  = res*stance  of  cell  heater  + lead  resistance  (50.060  ohms  + 0.286  ohm). 

The  current  source  was  monitored  over  a four-day  period  to  check  its  stability  and  accuracy.  A 
40-ohm  precision  resistance  box  and  a 10-ohm  NBS  standard  resistor  were  substituted  for  the  cell 
heaters.  Voltage  measurements  off  the  10-ohm  NBS  standard  resistor  were  made  with  the  ESI 
whenever  current  was  supplied  from  the  constant  current  source.  The  current  output  was  found 
to  be  0.8%  low  on  the  3 ma  range,  ±0.025%  to  ±0.29%  for  the  10  ma  range.  ±0.02%  to  ±0.05%  for 
the  60  ma  range  and  ±0.02%  to  ±0.05%  for  the  90  ma  range.  Current  stability  for  the  four-day  test 
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period  \v;is  ±0.33%  tor  the  3 ilia  range,  ±0.5%  tor  the  10  ma  range,  ±0.3%.  tor  the  60  ma  range  and 
*0.3'?  lor  the  90  ma  range.  The  aeeuraey  and  eurrent  stability  from  day  to  day  during  these  tests 
was  not  within  the  manufacturer's  specification  of  ±0.1%  and  ±0.05'/?  respectively.  This  presented 
no  problems  since  the  changes  made  to  the  calorimeter  allowed  current  measurements  for  each 
experiment. 

The  data  retrieval  system  was  augmented  with  the  addition  of  a high-frequency  filter 
(r  = 2s)  on  the  output  of  the  Keithley,  a Digitek  digital  voltmeter,  interface,  teletypewriter  and 
computer.  The  signal  from  the  digtal  voltmeter  is  fed  through  an  interface,  designed  by 
Mr.  Richard  Breen,  to  a computer  teletype  equipped  with  a paper  tape  punch  system.  The  paper 
tape  records  the  amplified  voltages  in  digitized  form  as  a function  of  time.  These  voltages  arc- 
processed  with  a computer  program  on  a Data-General-Nova-1 200  minicomputer  to  yield  the  area 
of  the  heat  curve.  This  method  of  data  retrieval  was  tried  because  occasional  errors,  in  the  heats 
measured,  as  much  as  ±2%  were  found  using  the  recorder’s  ball-and-disk  integrator.  The  program 
w hich  processes  the  voltage  data  of  an  experiment  can  be  found  in  appendix  A.  The  software  was 
tested  using  the  calibration  heater  system  of  the  calorimeter  to  provide  simulated  heats  of  dilution. 
Experiments  were  performed  to  determine  the  data  collection  rate  and  the  number  of  points 
necessary  for  an  accurate  heat  determination.  A data  collection  rate  of  one  voltage  point  every 
6 seconds  was  chosen  since  an  error  of  only  0.04%  of  the  experimental  heat  was  found  at  this 
setting  when  compared  with  a one-second-date-collection  rate  for  a 0.6-calorie  heat  effect.  The 
heat  curve  shown  in  figure  1 was  divided  into  three  sections.  The  points  in  the  first  and  third  section, 
K and  M respectively  were  used  to  calculate  a least  squares  baseline.  The  points  in  section  N were 
then  integrated  using  Simpson’s  rule  to  provide  the  curve  area  under  section  N.  The  area  between 
the  curve  and  baseline  is  then  calculated  as  the  area  under  the  curve  minus  the  area  under  the 
baseline.  The  number  cr  voltage  points  used  for  each  curve  section  was,  20  points  (2  min)  for  K. 

150  points  (15  min)  foi  N and  50  points  (5  min)  for  M.  These  values  were  chosen  since  they 
allowed  measurement  of  greater  than  99.9%  of  the  experimental  heat  and  provided  a long-term 
precision,  based  on  calibration  experiments  with  mixing  shown  in  tables  A-l  to  A-7,  appendix  A 
of:  ± 1 .8%  for  the  0.01  -niv  scale,  ± 1 .3%  for  the  0.03-mv  scale,  ±0.6%  for  the  0. 1 -mv  scale,  ±0.5% 
for  the  0.3-mv  scale,  ±0.2%  for  the  1 -mv  scale,  ±0. 1 % for  the  3-mv  scale  and  ±0.08%  for  the  1 0-mv 
scale.  Thus  the  software  retrieval  system  represents  an  order  of  magnitude  improvement  in  the 
precision  of  the  data  collected  with  the  ball-and-disk  integrator  of  the  recorder. 

The  results  of  tables  A-l  to  A-7,  appendix  A indicate  that  the  precision  of  calibration 
constants  determined  for  several  months  was  excellent  and  was  greater  on  high-heat  ranges. 

Because  of  these  results,  experiments  were  performed  to  determine  range  factors  which  would 
allow  calculation  of  calibration  constants  for  any  heat  range  from  a single  calibration  measurement 
on  the  10-mv  range  each  day.  The  range  factors  are  ratios  of  the  Keithley  output  voltage  to  the 
input  voltage.  These  were  measured  with  an  ESI  voltohnur.eter  off  an  NBS  standard  resistor  when 
the  Keithley  microvoltmeter  was  driven  full-scale.  A typical  circuit  diagram  for  these  experiments 
is  shown  in  figure  2.  The  results  of  two  ratio  experiments  performed  two  weeks  apart  were 
averaged  and  are  shown  in  table  A-8  (appendix  A)  as  range  factors.  Assuming  the  10-mv  range  of 
the  Keithley  to  be  correct  (range  factor  = 1 ),  the  following  adjusted-Keithley  ranges  (KR)  were 
obtained  by  dividing  the  respective  Keithley  range  by  its  range  factor:  the  3-mv  range  is  then 
2.9904  mv,  the  1-mv  range  is  0.9962  mv,  the  0.3-mv  range  is  0.3000  mv.  the  0.1 -mv  range  is 
0.1000  mv  and  the  0.03-mv  range  is  0.0301  mv.  Electrical  noise  problems  prevented  the  accurate 
measurement  of  the  0.01-mv  range  factor,  and  consequently  this  range  factor  was  set  at  unity, 
making  the  0.01-mv  Keithley  range  exactly  0.01  mv.  Based  on  the  other  range  measurements, 
setting  this  range  factor  at  unity  could  lead  to  calibration  constant  errors  for  this  range  of  ±0.4% 
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Figure  2.  Circuit  Diagram  Used  to  Determine  Keithley  Range  Factors 


which  is  acceptable  since  the  long-term  calibration  constant  precision  was  ± 1 .8%.  The  following 
equation  indicates  how  these  adjusted-Keithley  ranges  and  a single  10-mv  electrical  calibration  with 
mixing  each  day  were  used  to  calculate  calibration  constants  for  all  the  other  heat  ranges: 

ed  = vein X KR)/10mv  (8> 


where  e 

= calibration  constant  for  the  desired  Keithley  range 

e10mv  = 10_mv  calibration  constant  calculated  from  equation  7 

KR  = actual  Keithley  range  obtained  by  using  the  range  tactors. 

Table  A -8  presents  a comparison  between  the  average  calibration  constants  lor  various  heat 
ranges  obtained  from  long-term  calibration  experiments  listed  in  tables  A-l  through  A-7,  and 
calibration  constants  for  these  ranges  calculated  using  the  average  10-mv  calibration  constants  is 
good  since  the  largest  difference  noted  is  only  0.35%  for  the  1-mv  range.  This  procedure  was 
justified  by  the  fact  that  the  Keithley  range  resistors  are  rated  at  ±0.025'%  and  the  stability  of  the 
digital  voltmeter  and  ESI  are  ±0.01%  and  ±0.02  respectively.  Thus  only  one  calibration  per  day  on 
the  1 0-mv  range  provides  calibration  constants  for  all  the  heat  ranges  of  the  calorimeter.  1 his 
procedure  represents  a savings  in  time  ot  as  much  as  live  hours  per  day  depending  on  the  number  ol 
experiments  performed. 

On  a high-sensitivity  range  (1  pv),  electrical  noise,  errors  due  to  the  heat  of  friction  ot 
the  sample  liquid  against  the  cell  walls  and  baseline  drift  are  at  a maximum.  An  attempt  was  made 
to  improve  the  accuracy  of  the  heats  under  these  conditions  by  reducing  the  number  of  points 
measured.  Less  baseline  drift  would  be  recorded  and  any  loss  in  the  experimental  heat  measured 
could  be  corrected  with  a calibration  run  under  the  same  conditions.  Ten  electrical  calibrations 
with  mixing  were  performed  on  the  1-pv  range  by  inputing  1 ma  through  the  sample  cell  heater 
for  45  seconds.  The  areas  under  these  curves  were  calculated  using  several  different  curve-section 
times.  The  results  are  listed  in  table  A-9.  The  heat  curves  were  found  to  fall  into  two  classes  due  to 
different  heat  curve  profiles.  Six  of  the  electrical  calibrations  produced  normal  heat  curves  like  the 
one  shown  in  figure  3. A while  four  of  the  calibrations  had  a large  negative  peak  followed  by  a 
positive  peak  at  the  beginning  of  the  heat  curve  like  that  shown  in  figure  3.  B.  This  anomalous 
peak  in  four  of  the  ten  experiments  appears  to  be  a result  ot  mechanical  stress  on  the  semiconductor 
thermopiles,  and  the  area  averages  shown  for  these  curves  are  equivalent  to  the  effect  ot  a loss  ot 
approximately  30%  of  the  experimental  heat  and  show  a 1 5%  increase  in  the  standard  deviation  of 
a measurement.  It  should  be  pointed  out  that  the  lowest  heat  range  used  in  the  later  experimental 
work  was  10  pv  so  that  the  effect  is  one-tenth  as  large.  Whenever  it  occured  the  experiment  was 
rerun.  Based  on  the  average  curve  areas  for  the  six  calibrations  with  normal  heat  profiles  listed  in 
table  A-9.  the  optimal  curve  section  time  which  should  be  used  to  compute  curve  areas  on  the  1 -pv 
scale  is  2-10-4.  This  setting  detects  96%  of  the  experimental  heat  with  a standard  deviation  tor 
an  individual  measurement  of  2 volt  seconds  or  7 p calories.  Since  these  calibration  experiments 
were  performed  before  optimal  control  of  the  calorimeter  air  bath  was  achieved  and  without  the 
high-frequency  filter  on  the  output  of  the  Keithley  microvoltmeter,  another  series  of  calibrations 
were  performed  on  the  1-pv  scale  with  mixing  after  all  the  calorimeter  moditicationslhad  been 
finalized  in  order  to  measure  the  calorimeter’s  sensitivity.  The  areas  for  these  experiments  were 
computed  using  a curve  section  time  of  2-15-5  in  order  to  measure  100%  ot  the  heat.  The  areas 
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were  then  multiplied  by  a calibration  constant  calculated  from  the  average  10-mv  calibration 
constant  and  equation  8.  The  results,  presented  in  table  A-10,  indicate  that  the  sensitivity  of  the 
calorimeter,  given  by  the  standard  deviation  of  a single  heat  measurement  is  ±21  ^calories.  This 
value  agrees  exactly  with  that  obtained  when  the  standard  deviation  of  the  class  A area  for  a curve 
section  time  of  2-15-5  in  table  A-9  is  multiplied  by  the  same  calibration  constant  used  in  table  A- 1 0. 
This  value  also  agrees  favorably  with  the  value  of  ±20  ^calories,  reported  by  Falcone.39 

D.  Calorimetric  Procedure. 

The  solutions  were  made  up  on  a single-pan  Mettler  balance  capable  of  weighing  to 
±0. 1 mg.  The  sample  cell  was  rinsed  with  water,  followed  by  a rinse  with  Fischer-certified  methanol 
and  dried  with  a stream  of  nitrogen.  W(2)  grams  of  solution  containing  N(2)  moles  of  solute  at 
concentration  mj  (2),  and  W(4)  grams  of  solvent  water  were  loaded  by  weight  into  the  2-  and  4-cc 
cell  compartments  with  5-cc  plastic  syringes  equipped  with  lengths  of  Teflon  tubing.  When  a 
mixture  was  measured,  N(2)  is  the  moles  of  solute  a,  (Na),  plus  the  moles  of  solute  b,  (N^).  The 
materials  were  allowed  to  equilibrate  at  the  calorimeter  temperature  until  a stable  baseline  was 
achieved  on  a heat  range  one-hundred  times  more  sensitive  than  the  experimental  range  (1-2  hr). 

The  first  dilution  was  begun  by  rotating  the  calorimeter  unit  through  at  least  one  rotation  cycle, 
resulting  in  q calories  of  heat  being  released  and  a final  solution  concentration  of  mf.  A series  of 
rotations  were  performed  after  the  first  dilution  of  each  solute  until  there  was  no  detectable  heat 
effect.  The  number  of  rotations  for  each  solute  series  was  then  used  in  subsequent  dilutions  of 
these  solutes  to  assure  complete  mixing.  After  the  initial  baseline  was  reestablished,  a calibratiion 
(with  mixing)  was  performed  on  the  10-mv  range.  Keithley  range  factors  were  then  used  to  calculate 
calibration  constants  for  all  the  experimental  heat  ranges  from  the  10-mv  calibration.  Previously, 
calibrations  were  performed  after  every  dilution.  The  4-cc  compartment  was  then  emptied  by 
weight.  W(4)out  grams  was  withdrawn,  with  a 10-cc  plastic  syringe  leaving  16.5  ± 2 mg  of  solution 
on  the  walls  of  the  4-cc  cell  or  8.4  ± 2 mg  if  the  2-cc  compartment  was  emptied.  Pure  solvent 
W(4)jn  grams  was  then  loaded  into  the  4-cc  compartment  producing  a dilute  solution  containing 
N(4)  moles  of  solute  at  concentration  mj  (4).  This  dilute  solution  concentration  was  calculated 
from  the  number  of  moles  of  solute  left  on  the  walls  N(4)  and  the  total  weight  of  solvent  added  to 
the  4-cc  compartment.  Thus  the  ensuing  dilutions  involve  N(4)  moles  of  solute  at  concentration 
mj  (4)  mixing  with  N(2)  moles  of  solute  at  concentration  mj  (2)  to  give  N(2)  + N(4)  moles  of  solute 
at  concentration  mf'.  The  concentration  of  the  solution  in  either  cell  is  easily  calculated,  since 
the  weights  of  solution  and  solvent  added  and  withdrawn  and  the  amount  of  solution  retained  on 
the  walls  is  known  at  any  point  in  a series  of  dilutions.  The  following  equations  indicate  the 
computer  procedure  used  to  calculate  the  solute  concentration  in  either  cell  for  a series  of  dilutions. 

Initially  the  2-cc  cell  is  loaded  with  W(2)  grams  of  solution  at  molality  m,  (2)  and  the 
4-cc  cell  is  filled  with  W(4)  grams  of  water. 

WTA  = W(2)/[  1000.  + MW  X mj(2)] 

WTBB  = W(4)  X 0.001 


N(2)  = mj  (2)  X WTA 

mf  = N(2)/(WTA  + WTBB) 


where 

WTA  = the  number  of  grams  of  solvent  in  W(2)  grams  of  solution 
WTBB  = the  number  of  kg  of  solvent  used  to  dilute  W(2)  grams  of  solution 
MW  = molecular  weight  of  solute. 

At  this  stage,  the  contents  of  solution  in  cell  2 are  mixed  with  the  solvent  in  cell  4, 
creating  a uniform  solution  concentration  mf  in  each  of  the  two  cells.  The  4-cc  compartment  is 
then  emptied  [W(4)out  grams  being  removed]  leaving  16.5  ± 2 mg  of  solution  ofmf  on  the  cell 
walls.  W(4)jn  grams  of  solvent  are  then  added,  creating  a dilute  solution  mj(4)  in  the  4-cc  cell. 
Another  rotation  is  then  performed  resulting  in  a new  solution  concentration  mf ' in  both  cells. 

mj  (2)  = mf  (from  previous  dilution) 

[W(2)  + W(4)]  - [W(4)out  + 0.0165] 

WTA  = 

1000  + MW  X mf 

WTB  = 0.0165/(1.000.  + MW  X mf) 

WTBB  = [W(4)in]  X 0.001 

N(2)  = mf  X WTA 

N(4)  = mf  X WTB 

mi(4)  = N(4)/(WTB  + WTBB) 

mf ' = ]N(2)  + N(4)/WTA  + WTB  + WTBB) 

where 

WTB  = the  number  of  grams  of  solvent  remaining  in  the 
4-cc  cell  after  it  was  emptied. 

Dilutions  were  performed  in  this  manner  until  the  final  concentration  mf  was  0.1  molal  or  less. 

The  experimental  heat,  q,  is  related  to  the  initial  and  final  solution  excess  enthalpies  by 
the  following  equation: 


> 


where  Hex  (m)  is  the  excess  enthalpy  of  the  solution  per  mole  of  solute  at  molality  iu(</>l  - Hex). 
The  excess  enthalpy  of  the  solute  represented  by  the  equation: 

' Hex  = <t> l = Bjm  + B2Q}^  + B^m^  + . . . (10) 

was  substituted  into  equation  9 in  order  to  fit  the  experimentally  measured  values  for  the  change 
in  the  excess  enthalpy  per  mole  AHex  to  the  following  equation: 


-q/N(2)  - Hex  = j=(  n Bj 

where  N(F)  = total  moles  of  solute.  This  equation  was  used  in  a least  squares  procedure  to 
determine  the  coefficients  Bj  of  the  polynomial  and  to  extrapolate  the  heat  data  to  infinite  dilution 
The  right  hand  side  of  the  equation  is  the  heat  per  mole  of  solute.  Statistical  tests  were  performed 
to  assure  using  the  minimum  number  of  coefficients  to  represent  the  data.  Only  coefficients  which 
were  not  zero  with  95%  confidence  were  used. 

At  this  point  a brief  derivation  based  on  the  solution  theory  of  McMillan-Mayer  is  shown 
to  indicate  how  these  Bj  coefficients  were  used  to  calculate  the  pairwise  heats  of  interaction  | 
between  nonelectrolyte  solutes.  Wood1 4 has  shown  that  for  a system  containing  solutes  and  1 kg 
of  solvent,  the  excess  free  energy  Gex  may  be  expressed  as  a polynomial  in  the  molalities  of  the 
solute  species: 


‘N(F)m.j  - N(2)mj(2)>  - N(4)jiii(4)-»  * 


= AnlNiNjfg®iOJj  + AnlNiNjNk}g^imjfflk  + ... 


where  the  coefficient  A for  an  n particle  term  involving  nj  particles  of  species  Nj , n2  particles  of 
species  N2  etc.  is: 


n (nj !n2!n3!  . . .) 

These  An  coefficients  permit  {NjNj}  = {NjN:J  = {NjNj  J for  the  three  cases  in  which  species  i and  j 
are  the  same  except  for  some  type  of  label  which  does  not  affect  the  intermolecular  forces  between 
these  molecules.  The  species  within  the  brackets  {NjNj  | R correspond  to  the  McMillan-Mayer1 6 
cluster  integrals  and  indicates  that  the  cluster  integrals  are  functions  of  the  excess  free  energy.  The 
equation  for  the  excess  heat  content  is  obtained  by  taking  the  partial  derivative  of  the  excess  free 
energy  with  respect  to  the  reciprocal  absolute  temperature  as  follows: 


[NiNi]h  - 


6 ( [NjNj  ] g/T) 
61/T 


6 (Gw/D  = Hex 
61/T 
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By  applying  this  derivative  to  equation  (12)  a power  series  in  molality  results: 


Hex  , 

— - AnlNjNjlbEjjjij  + AnlNjNjNklbiDiinjnik  + .. 


where  the  subscript  h indicates  an  enthalpy  of  interaction.  Since  all  of  the  work  performed  in  this 
study  involved  the  enthalpy  of  dilution  of  single  and  binary  nonelectrolyte  solutions,  equation  1 3 
was  used  to  express  the  change  in  excess  enthalpy  per  mole  for  a dilution  of  a solution  containing 
Xsmj  molal  sucrose  (S)  and  xun}j  molal  urea  (U)  as  follows: 

~^T=  (mrmf) [2XSXU | US th+*u  I UU fh+xs1 SS } h1  + XjXu  (mf-m2)  [3{lJSS}h] 

+ XgX^WrWf)  [3|UUS}h]  + Xs(mj2-mf)[{sSS}h]  + Xu(mf-m2)[{uuu}h]  + ...  (14) 

where  n is  the  total  moles  of  solute  of  all  types,  mj  the  initial  molality,  mf  the  final  molality  after  a 
dilution  and  x the  mole  fraction  of  the  respective  solute.  Substitution  of  the  experimental  heat  data 
data  for  this  dilution  into  the  least  squares  equation  (11),  results  in  a Bj  coefficient  which  is  equal 
to  the  first  term  of  equation  (14):  2xsxu  {US  } h + x^  {UU}  h + xf  |SSj  h.  Thus  this  Bj 
coefficient  contains  all  of  the  pairwise  heats  of  interaction  which  result  when  a binary  urea-sucrose 
solution  is  diluted.  The  dilution  of  a single  nonelectrolyte  solution  results  in  a Bj  coefficient  which 
directly  gives  the  heat  of  interaction  between  the  two  like  solute  molecules  since  in  this  case,  using 
the  dilution  of  a urea  solution  as  an  example,  the  change  in  excess  enthalpy  per  mole  is  given  by: 

= {UUfh(mrmf)  + |UUU}h(m?-m2)  + (15) 

Thus  the  heat  of  interaction  between  to  unlike  nonelectrolyte  solute  molecules  was  calculated 
from  the  knowledge  of  the  like  pairwise  heats  of  interaction  obtained  by  diluting  single  nonelectrolyte 
solutions  and  combining  these  values  with  the  interaction  expression  obtained  for  the  dilution  of 
a binary  mixture  of  these  solutes.  Table  A-l  1,  Appendix  A gives  heats  of  dilution  values  in  cal/mol 
for  many  single  and  binary  nonelectroltye  solutes.  Table  A-l  2,  Appendix  A lists  the  sum  of  the  heats 
of  interaction  for  an  equal  mole  fraction,  binary  solution  of  urea  (U)  and  pentaerythritol  (P),  i.e., 
1/4(2  jUPfh  + {UUlb  + { PP { b)  as -10.5  cal  kg/mol^.  The  heats  of  interaction  for  the  single 
nonelectrolyte  solutions  are  also  reported  in  table  3 as  { UU  } ^ = -83.7  and  { PP  } h = 94.4  cal 
kg/mol*-.  Substituting  these  values  into  the  expression  for  the  1-1  urea  pentaerythritol  dilution 
gives -10.5  = l/4[2  { UP } b+  {-83.7}  + {94.4 } ] and  solving  for  { UP  | ^ yields  the  pairwise  heat  of 
interaction  between  a urea  and  pentaerythritol  molecule  { UP  } h = -26.4  cal  kg/mol^  as  shown  in 
table  A-l 3,  Appendix  A. 

An  approximately  equal  number  of  dilutions  were  performed  with  the  roles  of  the 
2-  and  4-cc  compartments  being  reversed.  The  experimental  dilutions  were  usually  continued  until 
the  heat  of  dilution  was  approximately  100  calories:  at  which  time,  the  sample  cells  were  emptied, 
rinsed,  dried  and  loaded  for  another  series  of  dilutions. 
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A correction  for  the  wet  wall  heat  of  friction  of  approximately  0 ± 20  /acal,  produced  by 
the  motion  of  the  liquid  sample  against  the  cell  walls  during  mixing,  was  made  whenever  the 
experimental  heat  of  dilution  was  60  meal  or  less.  The  method  used  for  this  correction  was  to 
subtract  from  the  heat  of  dilution,  the  heat  of  friction  observed  upon  rotation  of  the  homogeneous 
solution  produced  by  the  previous  dilution.  The  estimated  error  for  this  procedure  on  the  reported 
heat  values,  was  ±20  /a  calories. 

Since  the  heat  of  friction  of  an  initial  dilution  (cell  walls  dry)  has  been  reported  to  be  as 
much  as  ± 200  ^calories  by  Levine,3  8 a slightly  different  procedure  was  used  for  a few  dilution 
runs  which  gave  initial  heats  of  dilution  of  only  1 00  mealories  or  less.  For  these  runs,  the  4-cc  cell 
was  initially  filled  with  approximately  5 mg  of  solution  by  weight.  The  cells  were  then  rotated  and 
any  solution  found  in  the  2-cc  cell  was  withdrawn  by  weight  leaving  8.4  ± 2 mg  of  solution  on  the 
walls.  The  2-cc  cell  was  then  filled  with  a weighed  amount  of  solvent  and  the  normal  procedure 
for  the  dilutions  was  used  from  this  point.  The  dilution  runs  performed  by  this  procedure  are 
labelled  wet  wall  in  table  A-l  1 of  the  results  section  and  are  indicated  by  (C). 

All  of  the  single  and  binary  nonelectrolyte  dilutions  performed  in  this  work  were  plotted 
as  -AHex/(rrij-mf)n  versus  (mj+njf)  in  figures  4-10  of  the  results  section.  The  intercepts  of  these 
plots  gives  B,  of  equation  1 1 directly  as  the  sum  of  all  the  pairwise  heats  of  interaction,  and  the 
slope  of  these  curves  yields  a sum  containing  all  the  triplet  heats  of  interaction.  Since  we  were 
interested  in  only  the  interaction  heats  of  molecular  pairs,  it  was  important  to  measure  the  enthalpy 
of  dilution  at  low  concentrations.  At  concentrations  below  0.2  molal,  the  probability  of  molecules 
having  more  than  pair-type  interactions  becomes  quite  small.  Therefore  the  extrapolation  of 
-AHex/(mj  -mf)n  to  infinite  dilution  provides  the  pairwise  heats  of  interaction  with  good  accuracy. 

In  order  to  plot  the  data  in  the  form  -AHex/(mj-mf)n  versus  (mj+mf),  it  was  necessary 
to  correct  it  for  the  fact  that  the  solvent  compartment  (4-cc)  periodically  contained  small  amounts 
of  solute  in  addition  to  the  diluent  solvent.  Since  the  number  of  moles  of  solute  in  the  4-cc 
compartment  was  quite  small,  the  calculation  of  the  heat  that  would  have  been  observed  if  the  4-cc 
compartment  had  been  filled  with  pure  solvent  involved  only  a small  correction.  These  corrected 
AH®^  values  shown  in  table  A-l  1 of  the  results  section  for  the  nonelectrolyte  mixtures,  were 
obtained  by  using  the  experimental  concentration  data  for  the  2-  and  4-cc  compartments,  together 
with  the  experimental  heats  of  dilution  q,  in  the  following  equation: 


Mf^rr  = 0L(;orr  = * HeX[mi(2)]  = -q/N(2)  + N(4)/N(2){HeX[mi(4)]  -HeX(mf)}  (16) 


where  Hex[ijij(4)]  and  Hex(mf)  were  evaluated  from  the  coefficients  of  the  final  least  squares  fit 
of  the  data.  This  correction  never  exceeded  2%  of  q. 


III.  RESULTS. 


The  experimental  dilution  data  for  all  the  single  and  binary  nonelectrolyte  mixtures  were 
fitted  to  the  least  squares  equation  ( 1 1 ) previously  described.  A computer  program  used 
equation  (1 1 ) to  express  the  change  in  excess  enthalpy  per  mole  as  a function  of  the  initial  and 
final  molalities  of  the  solutions  being  diluted.  The  following  information  from  the  computer 
outputs  is  presented  in  table  A-l  1,  appendix  A: 


rrij(2),  N(2),  N(4),  mj- defined  earlier 
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KHKfc 


Xa  = Na/(Na+Nb) 


mole  fraction  of  solute  a in  a mixture  (a+b)  excluding  solvent. 


HSp=  -^(2) 


= the  experimental  change  in  the  excess  enthalpy  per  mole  of 
the  solute  in  the  more  concentrated  solution. 


^corr  = exPerimental  change  in  excess  enthalpy  per  mole  of  the 
solute  in  the  more  concentrated  solution  equal  to  -q/N(2) 
corrected  for  setting  N(4)  equal  to  zero  (see  equation  16). 

cx 

^calc  = t*ie  calculated  least  squares  value  for  Hcorr  above. 


The  B;  coefficients  required  to  fit  the  heat  data  to  the  least  squares  equation  (1 1 ) are 
shown  in  table  A-12  together  with  the  standard  deviations  for  the  coefficients. 


As  shown  earlier,  for  a dilution  of  a single  nonelectrolyte  solute  A,  the  first  coefficient 
obtained  from  a least  squares  analysis  of  the  heat  data  is  simply  the  heat  of  interaction  between 
two  A molecules,  i.e.,  B j = { AA  J b;  while  the  dilution  of  an  equimolal  mixture  of  A and  C yields 
a coefficient  Bj  = 1/4  [2  { AC  }h  + { AA } h + { CC  }h] . These  Bj  coefficients  were 
subsequently  used  to  calculate  (table  A-13)  pairwise  heats  of  interaction  | AB  } b for  all  the 
nonelectrolytes  studied. 

Finally , a series  of  - ^/(mj  -mf)n  versus  (mj+  mf ) plots  are  shown  in  figures  4- 1 0 are 

both  the  single  and  binary  nonelectrolyte  dilutions.  The  corrected  experimental  excess-enthalpy 
points  are  shown  in  relation  to  the  least  square  calculated  curves. 


(m(  + njj)  mol/kg 

Figure  4.  AH®£rr  Plotted  Versus  (rpj  + mf),  for  EG-S,  P-S,  G-S, 
P,  E-G  and  G Solutes 

The  line  gives  values  calculated  by  the  least  squares  fit. 
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Figure  10.  AHcorr  Plotted  Versus  (nij  + mf),  for  Four  Binary 
Nonelectrolytes  with  Pentaerythritol  as  One  Component 


IV.  DISCUSSION. 

A comparison  of  the  present  work  with  the  literature  indicates  that  very  little  work 
has  been  done  in  measuring  the  thermodynamic  properties  of  aqueous  nonelectrolyte  mixtures. 

The  available  data  on  single  nonelectrolyte  solutions  has  been  reviewed  by  Kozak,  Knight  and 
Kauzmann.1  0 The  pairwise  heat  of  interaction  for  sucrose  1 37.9  ± (1.4)  cal  kg/moU  measured 
here  agrees  within  experimental  error  with  the  values  reported  by  Lange40  of  1 33.6  and  Gucker 
and  Pickard4  1 of  134.6  cal  kg/mol“.  The  most  recent  and  most  accurate  literature  data  for  the 
heat  of  interaction  for  urea  was  reported  by  Hamilton  and  Stokes.4  2 Their  data  was  fit  by  the 
same  computer  program  used  in  treating  the  present  data  and  resulted  in  a value  of  -83.2  ± (1.2)  cal 
kg/mol  ~ which  is  in  excellent  agreement  with  the  present  result  of  -83.7  ± (2.7)  cal  kg/mol“.  Values 
for  urea  of -84.4  ± (2.9)  by  Cassel  and  Wood1  4 and  -85.9  by  Gucker  and  Pickard4  ' also  agree  well 
with  our  results.  Lange40  reported  a heat  of  interaction  for  glucose  of  82.1  cal  kg/mol“  which 
essentially  is  in  exact  agreement  with  the  value  of  82.0  ± (2.3)  found  in  this  study.  A heat  of 
interaction  for  ethylene  glycol  of  92.5  cal  kg/moI~  reported  by  Lange* 2 is  somewhat  outside  the 
estimated  experimental  error  with  the  value  of  86.5  ±(1.0)  found  in  this  study. 

Friedman  and  Krishnan*  1 have  measured  the  enthalpies  of  dilution  for  a series  of 
alcohols  and  binary  alcohol  mixtures.  Unfortunately,  the  accuracy  of  this  data  is  poor  because  these 
authors  did  not  have  a calorimeter  designed  for  dilution  experiments.  Another  contributing  factor 
to  the  poor  accuracy  obtained  by  Friedman  and  Krishnan  was  the  lack  of  an  accurate  extrapolation 
ot  -AHcx/(nij  -mf)n  versus  (mj+nif)  for  their  data.  The  value  of  (nij+rDf)  was  approximately  equal 
to  one,  which  is  too  large  to  obtain  an  accurate  extrapolation.  Thus  the  pairwise  interaction  heats 
obtained  from  this  extrapolation  probably  contain  contributions  from  triplet  terms.  The  more 
accurate  enthalpy  of  dilution  data  for  a series  of  alcohols  reported  by  Lange*  240  and  Franks*  are 
listed  in  table  A-14  of  this  section. 

It  must  be  emphasized  that  the  heats  of  interaction,!  fh’  reported  here  are  not  equal  to 
the  second  virial  coefficient  in  the  McMillian-Mayer  theory  of  solutions.  Friedman*  1 has  given  an 
equation  from  which  one  can  show  that  the  free  energy  interaction  parameter,  j L,  is  related  to  the 
McMillan-Mayer  second  virial  coefficient  B*  as  follows:  e 


(17) 


* Franks.  F.  Private  communication.  1975. 


BAA  = !AA!gVw  + RT0v(A)  - l/2jV^ 

where 

O 

vw  is  the  volume  of  water  in  liters  per  kilogram, 

0V  is  the  apparent  molal  volume  of  the  solute  A at  infinite  dilution 
j is  a constant  for  water  at  25°C  equal  to  1 . 1 X 1 O'* 
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The  heat  of  interaction  parameter  then  follows  from  equation  (17),  if  one  neglects  the  last  term 


6BAA/T  o <*RT2{AA}gVw 

= |AAfhVw  - RT20e(A)  (18) 

6 1/T  RT 

Thus  the  heat  term  shown  in  equation  (18)  is  directly  proportional  to  the  temperature  variation  of 
the  solute-solute  intermolecular  interaction.  By  analyzing  the  respective  terms  in  equation  (18) 
for  urea,  sucrose  and  a urea-sucrose  mixture,  an  estimation  of  the  magnitudes  of  these  terms  can 
be  obtained.  The  free  energy  interaction  parameters  reported  by  Ellerton  and  Dunlop44  and  the 
interaction  heats  and  partial  molal  expansibilities  measured  here  for  urea  and  sucrose  were  substituted 
into  equation  (18).  The  magnitude  of  each  term  is  shown  in  table  A-15  for  the  different  solutes. 

The  values  listed  in  table  A-15  show  that  the  heat  of  interaction  terms  are  the  most  important  in 
determining  the  value  of  the  temperature  variation  of  the  solute-solute  intermolecular  interaction 

Sli*  J'T 

-if! i . It  is  also  observed  that  the  partial  molal  expansibility  term  may  represent  as  much  as  40% 

5 1/T 

of  the  total,  whereas  the  free  energy  term  is  much  smaller  than  either  the  heat  or  expansibility  terms 
and  has  the  same  sign  as  the  heat  term.  It  should  therefore  be  possible  to  use  the  pairwise  heats  of 
interaction  to  look  for  trends  in  nonelectrolyte  group  interactions,  develop  sign  rules  for  these 
interactions  and  determine  whether  or  not  these  group  effects  are  additive.  Friedman1 1 has  already 
reported  a similar  study  for  monofunctional  alcohols  in  which  he  reported  a correlation  between 
the  values  of  the  heats  of  interaction  for  the  alcohols  and  the  length  of  their  alkyl  chains. 

Ratcliff  et  al. ,45'47  have  also  published  a number  of  papers  in  which  a group  solution  model 
was  used  to  predict  the  heats  of  mixing  alcohols  and  hydrocarbons. 

The  compounds  used  in  our  experiments  were  chosen  to  study  group  interactions 
between  amides,  hydrocarbons  and  polyhydroxy  compounds.  As  seen  from  table  A-13,  appendix  A 
the  amides  form  a series  of  increasing-hydrocarbon-chain  length  while  keeping  the  amide  contribution 
constant;  i.e.,  N-methylformamide  (NMF)  has  one  CH3  group,  N-methylacetamide  (NMA)  has 
two  CH3  groups,  N-methylpropionamide  (NMP)  has  a methylene  group  and  two  methyl  groups 
and  N-butylacetamide  (NBA)  has  three  methylene  groups  and  two  methyl  groups.  In  this  analysis 
no  distinction  was  made  between  a methylene  and  a methyl  group.  The  polyhydroxy  compounds 
were  selected  for  an  increase  in  the  number  of  hydroxy  groups  per  molecule.  Thus  ethylene 
glycol  (EG)  has  two  OH  groups,  pentaerythritol  (P)  has  four  OH  groups,  and  glucose  (G)  and 
sucrose  (S)  each  have  five  and  eight  OH  groups  respectively.  Urea  was  also  included  to  represent 
a pure  amide  group  and  has  arbitrarily  been  equated  to  one  and  one  half  amide  groups  in  this  work. 
Using  these  three  different  classes  of  molecules,  the  pairwise  group  interactions  between 
hydrocarbons,  polyhydroxy  compounds  and  amides  with  themselves  as  well  as  the  combinations 
between  these  groups  have  been  determined. 

The  pairwise  heats  of  interaction  for  the  series  of  amides  NMF-NMF,  NMA-NMA, 
NMP-NMP  and  NBA-NBA,  was  used  to  determine  the  hydrocarbon-hydrocarbon  group  interaction. 

In  this  amide  series,  the  amide  portion  of  the  molecule  remains  constant  as  the  length  of  the 
hydrocarbon  chain  increases.  The  heats  of  interaction  for  this  series  increases  in  value  thereby 
indicating  that  a hydrocarbon-hydrocarbon  group  interaction  is  positive.  This  result  agrees  with 


the  work  of  Friedman  and  Krishnan1 1 and  Franks*  on  aqueous  alcohols  where  they  found  the 
unlike  pairwise  enthalpies  of  interaction  to  be  positive  and  show  an  increasing  trend  as  a function 
of  the  length  of  the  alkyl  chain  on  an  alcohol.  A positive  heat  of  interaction  and  dependence  on 
the  length  of  the  alkyl  chain  is  also  expected  for  hydrocarbon  groups  in  an  aqueous  solution,  based 
on  the  free  energy  experiments  of  proteins  by  Nemethy  and  Scheraga.3  These  authors  maintain 
that  nonpolar  groups  preferentially  interact  with  each  other  rather  than  water  in  an  aqueous 
environment  by  hydrophobic  bonding,  resulting  in  an  entropy-dominated  process  which  can  be 
used  to  explain  the  large  positive  free  energies  of  solution  for  hydrocarbons  and  related  compounds 
in  water. 

Based  on  the  pairwise  heat  of  interaction  for  urea  with  itself,  amide-amide  group 
interactions  are  negative.  This  is  consistent  with  the  results  reported  by  Hamilton  and  Stokes,4  2 
and  Schellman5  for  urea  in  which  they  also  measured  an  endothermic  enthalpy  of  dilution.  The 
series  formed  by  urea  interacting  with  N-methylformamide  (NMF),  N-methylacetamide  (NMA), 
N-methylpropionamide  (NMP)  and  N-butylacetamide  (NBA)  provides  a measure  of  the 
amide-hydrocarbon  group  interactions.  The  amide-amide  contribution  should  remain  constant 
for  this  series,  and  since  the  hydrocarbon  chains  are  being  increased  on  one  of  the  solute  pair 
molecules,  an  amide  hydrocarbon  interaction  which  is  positive  is  observed.  Another  indication  that 
the  sign  of  the  amide-hydrocarbon  interaction  is  positive  is  provided  by  the  series  NMF  interacting  with 
NMA,  NMP  and  NBA.  Again  the  amide-amide  contribution  is  expected  to  remain  constant  for 
this  series  while  the  hydrocarbon  chain  length  on  one  solute  is  increasing.  If  the  contribution 
due  to  the  hydrocarbon-hydrocarbon  interaction  is  small  relative  to  the  amide-hydrocarbon 
interaction  then  the  trend  in  the  interaction  heats  observed  for  this  series  suggests  that  the  sign  of 
the  amide-hydrocarbon  interaction  is  positive. 

Focusing  attention  on  the  polyhydroxy  compounds  shown  in  table  A-13,  appendix  A, 
specifically  the  series  of  pairs  ethylene  glycol  (EG),  pentaerythritol  (P),  glucose  (G)  and  sucrose  (S). 
shows  that  the  interaction  heats  are  approximately  constant  for  EG,  P and  G,  whereas  sucrose  (S) 
shows  a moderate  increase.  This  indicates  that  the  heats  of  interaction  for  hydroxy  groups  are 
small  and  probably  positive.  The  fact  that  the  pairwise  interaction  heat  for  sucrose  is  larger  than 
the  rest  of  this  series  could  be  a result  of  steric  effects  or  the  availability  of  the  hydroxy  groups  to 
interact.  The  similarity  between  the  heats  of  interaction  for  ethylene  glycol  (EG),  pentaerythritol  (P) 
and  glucose  (G)  pairs  is  maintained  when  these  molecules  interact  with  each  other.  The  values  for 
EG-P,  EG-G  and  P-G  are  approximately  the  same.  Sucrose  again  remains  unique  for  the  polyhydroxy 
compounds  when  it  interacts  with  the  other  compounds  in  this  series.  Sucrose  (S)  interacting  with 
ethylene  glycol  (EG),  pentaerythritol  (P)  and  glucose  (G)  forms  a fairly  regular  trend  in  the  heats, 
with  S-EG  having  the  largest  pairwise  interaction  heat,  (larger  than  S-S  pairs),  and  S-G  the  smallest 
value.  Thus  it  appears  that  the  OH-OH  interaction  is  positive  and  probably  small  because  of  the 
equality  of  the  interaction  heats  between  molecules  with  two.  four  and  five  hydroxy  groups.  The 
fact  that  a fairly  regular  trend  exists  for  the  interaction  of  the  polyhydroxy  compounds  with 
sucrose  and  that  EG-S  is  the  largest  pairwise  interaction  heat,  suggests  a steric  mechanism  controlling 
the  ability  of  the  hydroxy  groups  to  interact  with  each  other. 

The  sign  of  the  hydrocarbon-hydroxy  group  interaction  heat  is  positive  based  on  the 
heats  of  interaction  for  the  series  of  amides  with  each  polyhydroxy  compound.  As  the  length  of 
the  amide  chain  is  increased,  the  heats  of  interaction  for  the  amide-polyhydroxy  pair  increase. 

This  series  also  shows  much  more  of  a regular  increase  in  the  interaction  heats  than  the  polyhydroxy 
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compounds  with  themselves.  As  one  increases  either  the  hydrocarbon  chain  or  the  number  of 
hydroxy  groups  available  for  interaction,  the  pairwise  heats  of  interaction  increase,  although 
pentaerythritol  seems  to  give  lower  values  than  expected. 

The  final  group  interaction  obtained  from  the  data  in  table  A-13.  appendix  A is  an 
amide  group  with  a hydroxy  group.  The  series,  urea  interacting  with  the  polyhydroxy  compounds, 
indicates  that  this  group  interaction  is  negative  since  urea-ethylene  glycol  is  slightly  positive  and 
urea  interacting  with  progressively  higher  hydroxy  substituted  molecules  becomes  increasingly 
more  negative.  It  is  disturbing  to  note  that  N-methylformamide  (NMF)  has  a pairwise  heat  of 
interaction  ot  65  cal  kg/mol~  while  N-methylacetamide  (NMA)  pairs  gjve  56  cal  kg/mol*-,  yet 
{ NMF-NMA } h = 88  cal  kg/niol^.  Based  on  simple  group  additivity,  { NMF-NMA  } b is  expected 
to  be  an  average  of  the  like  pair  interactions,  and  this  is  clearly  not  the  case.  Also,  the  interaction 
value  of  NMA  with  N-methylpropionamide  (NMP)  or  N-butylacetamide  (NBA)  is  consistently 
smaller  than  the  interaction  values  of  NMF  with  these  compounds.  This  is  in  spite  of  the  fact  that 
the  NMA  has  one  more  methyl  group  than  NMF  with  which  to  interact. 

Although  the  data  listed  in  table  A-l  3,  appendix  A contain  nonadditivities,  it  appears 
that  at  least  the  signs  of  the  various  group  interactions  have  been  obtained.  These  sign  rules  are 
shown  in  table  A-l  6,  appendix  A.  It  had  been  hoped  that  group  additivity  in  the  heats  of 
interaction  of  nonelectrolytes  would  be  found  in  this  work  since  the  work  of  other  authors  suggests 
this  to  be  the  case.  Friedman  and  Krishnan1 1 and  Franks,*  for  example,  have  shown  a correlation 
between  the  length  of  the  alkyl  chain  for  simple  alcohols  and  their  heats  of  interaction.  Schrier48 
has  also  shown  a group  additivity  for  amide  groups,  methyl  groups,  and  methylene  groups  interacting 
with  ions  from  the  salting-out  behavior  of  various  amides.  In  light  of  their  results,  the  heats  of 
interaction  shown  in  table  A-13,  appendix  A.  were  fitted  to  a linear  least  square  equation  to 
determine  if  a simple  additivity  existed  in  the  pairwise  heats  of  interaction  and  whether  a knowledge 
of  the  additive  interactions  could  be  used  to  predict  the  heats  of  interaction  of  nonelectrolytes 
which  have  not  been  measured.  We  assumed  that  every  group  on  a solute  molecule  could  freely 
interact  without  steric  hindrance  with  every  group  on  another  solute  molecule,  and  that  the 
interaction  between  groups  was  independent  of  nearest  neighbor  effects.  The  following  general 
equation  was  used: 


Wh  = C0  + (nCH2a  * nCH2b)ECH2-CH2  + (nAa  X nAjEA_A  ♦ 

("CHOHg  x nCHOHb)ECHOH-CHOH  + ("CH2a  x "Ab  + n(H2b  x nAa)ECH2  A 
+ (nCH2a  X nCHOHb  + "CH2b  x nrHOHa)tCH2-THOH  ♦ 

(nCHOHa  x nA„  * "CHOHb  x "aJ^HOHa  (19) 

where  nCH->a 's  nurT,ber  of  hydrocarbon  groups  on  a molecule  a.  nA^  is  the  number  of  amide 
groups  on  molecule  b.  is  the  number  of  poly  hydroxy  groups  on  molecule  a and  F^  ^ 

* Franks.  F.  Private  communication.  1975. 
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Ea.a.  etc-  are  interaction  energies  between  each  group  given  in  cal  kg/moK  The  term  Cq 
represents  the  interaction  between  a terminal  hydrogen  atom  (group)  on  a particular  solute 
molecule  with  the  available  hydrocarbon,  amide  and  polyhydroxy  groups  on  another  solute 
molecule.  If  each  of  these  interactions  were  included  in  equation  19,  Cq  would  have  been 
expressed  as  follows: 

C0  = (nHa  x nHb)EH-H  + (nHa  X nCH25  + nHb  X nCH2a)EH-CH2  + 

(nHa  X "Ab  + nHb  X nAa)EH-A  + (nHa  X nCHOHb  + nHb  X nCHOHa  )EH-CHOH 


The  interaction  of  a H atom  with  the  other  available  groups  was  not  expected  to  be  as  important 
as  the  interactions  between  the  groups  previously  defined  because  of  the  size  disparity  between  a 
H atom  and  these  groups.  The  inclusion  of  these  terms  would  also  have  resulted  in  more  variables 
than  were  warranted  by  the  data.  These  interactions  were  consequently  summed  in  the  form  of  a 
constant  Cq.  It  should  also  be  observed  that  the  hydroxy  group  previously  used  in  determining 
the  sign  of  group  interactions  has  been  expanded  to  include  the  methylene  carbon  and  hdyrogen 
as  part  of  this  group.  This  is  done  for  convenience  in  treating  the  heats  of  interaction  of  the  sugars. 

As  a result  of  this  analysis,  the  following  group  interaction  values  (E)  in  cal  kg/mol^ 
were  determined: 

ECH2-CH2  = 7-6  ± 5,  EA.A  = -88  ± 32,  ECHOH<:HoH  = 0.3  ± 0.5,  E<:H2-CHOH  = 9 1 1 • 
ECH2-A  = 14  ± 1 1 and  ECH0H.A  = -1 2 ± 4. 

The  plus  or  minus  limits  shown  for  each  group  interaction  value  represents  the  95%  confidence 
limits  for  these  results.  The  signs  of  the  interaction  heats  obtained  with  equation  19  agree  with 
our  earlier  observations  shown  in  table  A-l  6,  appendix  A.  Even  the  magnitude  of  the  CHOH-CHOH 
interaction  agrees  with  our  earlieT  estimates.  The  pairwise  heats  of  interaction  calculated  by  using 
these  group  interaction  parameters  are  compared  with  the  experimentally  determined  interaction 
heats  for  all  the  compounds  studied  in  this  work  in  table  A-l  7,  appendix  A.  An  indication  that 
these  analytical  interaction  values  could  only  be  obtained  by  a comprehensive  study,  such  as  that 
done  here,  is  seen  from  the  inability  to  obtain  statistically  significant  values  when  the  amide  series 
or  the  polyhydroxy  series  alone  was  analyzed  by  the  least-squares  equation. 

Although  the  standard  deviation  of  the  predicted  pairwise  interaction  heats  for  the 
compounds  shown  in  table  A-l  7 was  approximately  ± 45  cal  kg/mol“.  the  range  of  the  heats  of 
interaction  measured  extends  from  -100  to  over  +400  cal  kg/mol^  as  illustrated  in  figure  1 1 . This 
figure  illustrates  the  pairwise  interaction  heats,  calculated  using  the  group  interaction  parameters  as 
points,  versus  the  experimentally  determined  heats  of  interaction,  shown  as  a line  for  all  the 
nonelectrolytes  used  in  this  work,  as  well  as  several  other  compounds  reported  in  the  literature. 
These  results  provide  the  first  group  interaction  values  for  aqueous  binary  nonelectrolytes  at  low 
concentrations. 
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EXPERIMENTAL  {ABjh  cal  kg/mol2 

Figure  1 1 . The  Experimental  Versus  Calculated  Pairwise  Heats  of  Interaction 

• To  determine  if  the  group  interaction  values  could  be  used  to  calculate  the  pairwise 

heats  of  interaction  of  other  molecules  containing  amide  groups,  alkyl  groups  and  polyhydroxy 
groups,  the  enthalpies  of  interaction  for  the  alcohols,  reported  by  Franks*  and  Lange1  2 and  several 
other  nonelectrolytes,  were  compared  to  those  calculated  by  using  the  group  interaction  parameters. 

The  results  are  shown  in  table  A- 14.  The  agreement  between  the  calculated  and  literature  values  is 
as  good  as  the  least-squares  fit  of  our  data.  This  is  remarkable  in  light  of  the  number  of  simplifying 
assumptions  which  had  to  be  made  to  obtain  the  group  interaction  parameters.  The  effect  of 
structure  in  the  form  of  steric  effects  was  not  accounted  for  by  the  simple  group  additivity 
postulated.  Each  functional  group  on  one  solute  molecule  was  assumed  to  interact  with  every 
other  group  on  another  solute  molecule  no  matter  how  inaccessible  certain  groups  were  to  each 
other.  Because  of  the  inability  to  account  for  structural  effects  in  determining  various  group 
interaction  parameters,  we  would  expect  discrepancies  between  the  calculated  and  experimental 
heats  of  interaction  to  increase  as  the  amount  of  steric  hindrance  present  in  the  molecule  increases. 

• Franks,  F.  Private  communication.  1975. 

l . -2L. id 


This  accounts  for  the  identical  calculated  heats  of  interaction  obtained  for  n-butyl  alcohol  and 
t-butyl  alcohol  even  though  Franks*  has  reported  that  the  heat  of  interaction  of  t-butyl  alcohol  is 
lower  than  n-butyl  alcohol  because  of  steric  hindrance. 
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In  the  case  of  sugars,  both  glucose  and  sucrose  were  assumed  to  consist  bf  six  and  twelve 
polyhydroxy  groups  respectively.  One  would  expect  the  contribution  to  the  heat  of  interaction 
due  to  the  hemiacetal  and  hemiketal  oxygens  to  be  different  than  the  contribution  of  the  CHOH 
groups  on  the  molecules,  since  the  oxygen  atoms  are  totally  ring-bound,  whereas  the  CHOH  groups 
are  more  available  to  interact.  In  view  of  these  simplifying  assumptions,  the  calculated  results 
presented  in  table  A-14  are  quite  encouraging.  Except  for  isopentanol,  the  pairwise  heats  of 
interaction  of  the  remaining  compounds  have  been  calculated  to  within  45  cal  kg/moF. 

V.  CONCLUSIONS. 

This  work  was  originally  performed  to  determine  if  structure  rules  could  be  formulated 
and  used  to  predict  the  sign  of  AH  when  nonelectrolytes  and  their  mixtures  were  diluted  with 
water.  The  enthalpies  of  dilution  of  four  amides  (N-methylformamide,  N-methylacetamide, 
N-methylpropionamide  and  N-butylacetamide),  four  polyhydroxy  compounds  (ethylene  glycol, 
pentaerythritol,  glucose  and  sucrose),  and  of  urea  were  measured  in  this  work.  These  nonelectrolytes 
and  their  binary  mixtures  were  investigated  at  25°C  using  a modified  LKB  batch  microcalorimeter. 
The  enthalpy  data  for  these  solutes  was  used  to  study  the  pairwise  heats  of  interaction  between 
hydrocarbon,  amide  and  CHOH  groups.  Trends  found  in  the  pairwise  heats  of  interaction  data  for 
these  nonelectrolyte  solutions  resulted  in  the  following  sign  rules  for  the  group  heats  of  interaction: 

hydrocarbon-hydrocarbon  = +,  amide-amide  = -,  CHOH-CHOH  = + 
hydrocarbon-amide  = +,  hydrocarbon-CHOH  = +,  and  amide-CHOH  = -. 

The  concept  of  group  additivity  in  the  heats  of  interaction  of  nonelectrolytes  was 
' employed  to  fit  the  pairwise  heats  of  interaction  measured  in  this  work  to  the  following  linear 

least-squares  equation: 

|ab}h  = CQ  + £ n,Anf  Ejj- 
i,j 

It  was  assumed  that  every  group  on  a solute  molecule  could  freely  interact  without  steric  hindrance 
with  every  group  on  another  solute  molecule  and  a'so  that  the  interaction  between  groups  was 
independent  of  nearest  neighbor  effects.  For  the  first  time,  as  a result  of  this  work,  quantitative 
values  have  been  obtained  for  the  heats  of  interaction  between  hydrocarbon,  amide  and  CHOH 
groups.  The  standard  deviation  of  the  pairwise  interaction  heats  calculated  using  these  group 
interaction  parameters  was  approximately  45  cal  kg/mol2.  The  validity  of  these  results  was 
demonstrated  by  using  the  group  interaction  parameters  to  successfully  calculate  the  pairwise  heats 
of  interaction  of  the  following  previously  measured  aqueous  solutions:  methanol,  ethanol,  n-propanol, 
n-butanol,  t-butanol,  isopentanol,  glycerol,  glycolamide,  ethanol-n-propanol,  ethanol-n-butanol. 
ethanol-t-butanol,  n-propanol-n-butanol,  and  n-propanol-t -butanol. 


* Franks.  F.  Private  communication.  1975. 
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The  magnitude  of  the  group  heats  of  interaction  obtained  were  larger  than  the  magnitude 
of  the  hydrocarbon-hydrocarbon  heats,  except  for  the  CHOH-CHOH  interaction.  This  result 
is  significant,  since  the  hydrophobic  effect  has  been  used  to  explain  the  structure  of  protein  in  aqueous 
solutions.2  3 Yet  in  terms  of  the  group  heats  of  interaction,  amide-amide,  amide-hydrocarbon, 
amide-CHOH  and  hydrocarbon-CHOH  interactions  have  been  shown  to  be  more  important  than 
hydrocarbon-hydrocarbon  interactions. 

These  group  interaction  parameters  should  be  used  to  evaluate  the  pairwise  heats  of 
interaction  of  nonelectrolyte  systems  which  have  already  been  measured  and  predict  those  for  which 
there  are  no  data.  Work  should  also  be  performed  to  evaluate  the  assumption  that  interactions 
between  groups  are  independent  of  nearest  neighbor  effects.  A good  compound  to  study  would  be 
N-methylbutyramide  since  the  enthalpy  of  dilution  of  N-butylacetamide  has  already  been  measured 
in  this  work.  A comparison  of  the  results  for  these  solutes  should  determine  the  importance  of  the 
position  of  the  alkyl  groups  in  an  amide  molecule  (carbonyl  or  nitrogen  substituted)  on  the 
hydrocarbon-amide  interaction.  Additional  work  is  also  needed  to  more  accurately  assign  the 
magnitude  of  the  amide-amide  group  interaction.  The  enthalpy  of  dilution  of  formamide  should 
provide  this  information. 

In  order  to  investigate  the  thermodynamics  of  biological  systems,  the  heats  of  interaction 
for  amine  and  carboxylic  acid  groups  should  be  studied  in  acidic,  neutral  and  basic  pH  solutions.  If 
additivity  in  the  heats  of  interaction  for  these  groups  can  be  established,  it  should  be  possible  to 
predict  the  pairwise  heats  of  interaction  of  complex  biological  molecules  that  are  composed  of 
pieces  for  which  the  interactions  are  known. 
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APPENDIX  A 
TABLES 


Table  A-l . Long-Term  Calibration  Constants,  10  mv  Range  (with  Rotation) 


Date 

1/ma 

A/vs* 

e/mcal/vs 

3/28 

89.7966 

156.84 

36.908 

3/28 

89.7986 

156.67 

36.949 

4/2 

89.8125 

156.77 

36.938 

4/2 

89.8085 

156.75 

36.938 

4/4 

89.8502 

156.89 

36.940 

4/8 

89.8721 

157.07 

36.915 

4/8 

89.8741 

156.78 

36.985 

4/10 

89.8761 

156.95 

36.947 

4/10 

89.8780 

157.06 

36.923 

4/11 

89.8761 

156.89 

36.961 

4/11 

89.8800 

156.75 

- 36.998 

5/22 

89.9595 

157.06 

36.990 

5/22 

89.9992 

157.39 

36.945 

5/23 

90.0032 

157.16 

37.001 

5/24 

90.0131 

157.28 

36.981 

5/30 

90.0091 

157.40 

36.951 

5/30 

90.0032 

157.39 

36.948 

5/30 

90.0052 

157.35 

36.960 

5/30 

90.0072 

157.38 

36.952 

5/30 

90.0091 

157.39 

36.953 

5/31 

89.9923 

157.31 

36.958 

6/3 

90.0131 

157.03 

37.042 

6/4 

90.0111 

157.32 

36.970 

6/6 

90.0191 

157.46 

36.946 

6/7 

90.0210 

157.16 

37.017 

6/12 

90.0290 

157.46 

36.951 

6/14 

90.0330 

157.36 

36.980 

6/19 

90.0330 

157.60 

36.924 

6/21 

90.0290 

157.44 

36.958 

6/24 

90.0389 

157.40 

36.975 
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Table  A- 1.  (Contd) 


Date 

I/ma 

A/ vs* 

e/mcal/vs 

6/28 

90.0350 

157.40 

36.972 

. 7/1 

90.0389 

157.25 

37.01 1 

7/2 

90.0350 

157.59 

36.928 

' 7/3 

90.0310 

157.39 

36.971 

7/10 

90.0369 

157.29 

36.999 

7/10 

90.0369 

157.53 

36.942 

7/10 

90.0369 

157.39 

36.976 

7/10 

90.0369 

157.23 

37.013 

7/10 

90  0369 

157.31 

36.994 

7/10 

90.0369 

157.19 

37.023 

7/15 

90.0369 

157.26 

37.005 

7/16 

90.0429 

157.36 

36.988 

7/18 

90.0429 

157.33 

36.995 

7/19 

90.0389 

157.39 

36.977 

7/22 

90.0449 

157.49 

36.959 

7/23 

90.0429 

157.39 

36.980 

7/26 

90.0429 

157.24 

37.017 

7/31 

90.0429 

157.51 

36.953 

8/1 

90.0449 

157.35 

36.992 

8/2 

90.0409 

157.43 

36.970 

8/5 

90.0429 

157.40 

36.979 

8/6 

90.0449 

157.52 

36.951 

8/7 

90.0469 

157.43 

36.975 

8/8 

90.0449 

157.44 

36.970 

8/12 

90.0429 

157.47 

36.962 

8/13 

90.0449 

157.39 

36.983 

8/14 

90.0449 

157.25 

37.015 

8/30 

90.0072 

157.21 

36.994 

9/3 

90.0270 

157.25 

37.000 

9/4 

90.0330 

157.47 

36.954 

9/6 

90.0330 

157.27 

37.001 

Average  e = 36.971  mcal/vs 


a (Single  measurement)  = 0.029  = 0.08% 


* In  each  experiment  the  current  was  passed  for  exactly  60  seconds. 
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Table  A-2.  Long-Term  Calibration  Constants,  3 mv  Range  (with  Rotation) 


Date 

t/s 

I/ma 

A/ vs 

e/mcal/vs 

10/25 

11 

88.7062 

93.70 

1 1 .097 

10/26 

13 

89.2325 

111.79 

11.079 

10/29 

14 

88.7678 

119.07 

1 1.085 

12/14 

22 

89.5702 

190.64 

11.077 

12/14 

14 

89.5741 

121.26 

1 1 .084 

12/15 

18 

89.6059 

156.13 

11.075 

12/15 

18 

88.8571 

153.26 

11.095 

12/19 

18 

89.5165 

155.68 

11.085 

12/20 

18 

88.5671 

152.45 

11.081 

12/28 

18 

89.3080 

154.98 

11.083 

2/20 

18 

89.8483 

156.78 

11.090 

2/20 

20 

89.8105 

173.99 

1 1 .094 

2/20 

20 

89.8522 

174.19 

11.091 

2/21 

18 

89.8681 

156.86 

11.088 

2/22 

18 

89.9297 

157.17 

11.082 

2/29 

18 

89.8721 

156.44 

11.119 

Average  e = 1 1 .088  mcal/vs 

0 

(Single  measurement) 

= 0.010  = 0.1^ 
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Table  A-3.  (Contd) 


Appendix  A 


Date 

t/s 

l/ma 

A/vs 

e/mcal/vs 

1/2 

20 

49.7080 

160.13 

3.692 

1/2 

20 

49.7080 

160.12 

3.692 

1/2 

20 

49.7100 

160.31 

3.689 

1/2 

20 

49.7100 

160.31 

3.688 

1/2 

20 

49.7100 

160.44 

3.686 

1/17 

24 

49.7398 

191.67 

3.707 

1/17 

26 

49.7398 

208.45 

3.692 

1/17 

20 

49.7398 

160.34 

3.692 

1/21 

25 

49.7477 

200.66 

3.689 

1/21 

20 

49.7497 

160.45 

3.691 

1/21 

24 

49.7616 

192.68 

3.690 

1/30 

20 

49.7001 

16d.06 

3.693 

1/30 

20 

49.0883 

155.40 

3.710 

1/31 

20 

49.1876 

156.66 

3.696 

1/31 

20 

49.0605 

155.96 

3.693 

2/5 

20 

49.9364 

161.77 

3.689 

2/5 

20 

49.9364 

161.60 

3.692 

2/7 

20 

49.9206 

160.66 

3.712 

2/7 

23 

49.9265 

184.89 

3.710 

2/8 

20 

49.8530 

160.26 

3.711 

2/12 

20 

49.9345 

160.83 

3.710 

2/13 

20 

49.9285 

160.55 

3.716 

2/13 

20 

49.9285 

160.96 

3.706 

2/14 

20 

49.9145 

160.45 

3.716 

2/14 

20 

49.9225 

161.26 

3.698 

2/14 

20 

49.9265 

161.36 

3.697 

2/15 

20 

49.9305 

161.24 

3.700 

2/15 

20 

49.9305 

161.33 

3.698 

2/19 

20 

49.9305 

161.52 

3.694 

2/21 

20 

49.9682 

161.60 

3.697 

2/21 

20 

49.9742 

161.33 

3.704 

2/22 

20 

49.9702 

161.40 

3.702 

2/22 

20 

49.9642 

161.55 

3.698 

2/25 

20 

49.9881 

162.20 

3.686 

2/25 

20 

49.9881 

161.82 

3.695 

2/26 

20 

49.9901 

161.96 

3.692 

2/26 

20 

49.9901 

161.75 

3.697 

2/28 

20 

49.9901 

162.05 

3.690 

Average  e = 3.696  mcal/vs 
o (Single  measurement  = 8 X 10"^  = 0. 


49 


Table  A-4.  Long-Term  Calibration  Constants,  0.3  mv  Range  (with  Rotation) 


Date 

t/s 

1/ma 

A/vs 

e/mcal/vs 

10/5 

17 

— 

29.9587 

165.9 

1.101 

10/5 

20 

29.9586 

194.6 

1.104 

10/18 

17 

29.6727 

162.3 

1.103 

10/26 

12 

29.7879 

115.1 

1.106 

10/29 

15 

29.6449 

141.6 

1.114 

11/1 

16 

29.4403 

148.9 

1.1 14 

11/20 

21 

29.5158 

195.0 

1.123 

11/26 

19 

29.7025 

180.5 

1.111 

12/3 

20 

29.2079 

183.6 

1.112 

12/6 

22 

29.8395 

210.6 

1.113 

12/6 

13 

29.8395 

124.7 

1.111 

12/7 

21 

29.8415 

201.0 

1.113 

12/14 

17 

29.8730 

163.3 

1.111 

12/14 

21 

29.8772 

202.0 

1.1 10 

12/15 

17 

29.8812 

163.0 

1.114 

12/15 

17 

29.6270 

159.6 

1.118 

12/19 

18 

29.5158 

168.4 

1.1 14 

12/20 

17 

29.5714 

159.9 

1.112 

12/27 

20 

29.7183 

189.8 

1.113 

1/21 

17 

29.8593 

163.3 

1.111 

1/30 

17 

29.4423 

161.2 

1.094 

1/31 

17 

29.5078 

159.2 

1.112 

1/31 

20 

29.5256 

188.2 

1.108 

2/7 

17 

29.9587 

163.4 

1.118 

2/7 

17 

29.9130 

163.1 

1.116 

2/8 

17 

29.9388 

163.7 

1.114 

2/13 

19 

29.9627 

182.4 

1.119 

2/13 

17 

29.9567 

163.8 

1.114 

2/14 

18 

29.9607 

173.5 

1.114 

2/15 

19 

29.9646 

183.7 

1.112 

2/19 

19 

29.9627 

184.1 

1.108 

2/20 

17 

29.9587 

163.4 

1.118 

2/20 

19 

29.9527 

183.0 

1.1 14 

2/20 

17 

29.9646 

164.0 

1.113 

2/21 

20 

29.9646 

193.7 

1.109 

2/21 

17 

29.9906 

164.1 

1.115 

2/21 

17 

29.9924 

163.9 

1.1 16 

2/22 

17 

29.9880 

164.4 

1.112 

2/22 

17 

29.9905 

164.3 

1.113 

2/25 

17 

30.0004 

165.3 

1.108 

2/26 

17 

30.0004 

164.5 

1.112 

2/28 

17 

30.0004 

164.7 

1.111 

Average  e = 1.112  mcal/vs 

a ( 

Jingle  measurement)  = 

5 X 1 0‘3  = 0.5% 
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Table  A-5.  I.ong-Term  Calibration  Constants,  0.1  mv  Range  (with  Rotation) 


Date 

t/s 

1/ma 

A/ vs 

e/mcal/vs 

10/5 

16 

19.9817 

211 

0.362 

10/5 

9 

19.9817 

117 

0.366 

10/5 

12 

19.9797 

166 

0.371 

10/18 

16 

19.7910 

204 

0.367 

10/29 

14 

19.7732 

117 

0.370 

10/31 

11 

19.4300 

135 

0.369 

10/31 

11 

19.7811 

139 

0.371 

11/19 

14 

19.6540 

175 

0.370 

11/21 

18 

19.6361 

226 

0.368 

11/26 

18 

19.8070 

228 

0.370 

11/26 

12 

19.8030 

152 

0.371 

12/3 

17 

19.5150 

209 

0.371 

12/5 

11 

19.8308 

140 

0.369 

12/6 

11 

19.8983 

141 

0.370 

12/7 

18 

19.9003 

232 

0.368 

12/7 

10 

19.9003 

128 

0.369 

12/7 

12 

19.9023 

154 

0.368 

12/15 

17 

19.7632 

211 

0.376 

12/20 

14 

19.6262 

175 

0.370 

12/27 

15 

19.8228 

190 

0.370 

12/28 

17 

19.8665 

217 

0.370 

1/17 

18 

19.9062 

230 

0.370 

1/21 

14 

19.9122 

179 

0.371 

1/21 

17 

19.9142 

218 

0.370 

1/21 

13 

19.9102 

166 

0.372 

1/30 

13 

19.8904 

166 

0.371 

1/31 

17 

19.6897 

213 

0.370 

1/31 

13 

19.6898 

163 

0.369 

1/31 

14 

19.6262 

174 

0.370 

2/5 

15 

19.9837 

194 

0.369 

2/5 

17 

19.9837 

218 

0.372 

2/8 

13 

19.9678 

166 

0.373 

2/12 

13 

19.9837 

168 

0.371 

2/14 

13 

19.9817 

167 

0.372 

2/14 

16 

19.9797 

206 

0.371 

2/14 

13 

19.9797 

169 

0.368 

2/15 

17 

19.9817 

219 

0.370 

2/15 

13 

19.9817 

167 

0.371 

2/19 

17 

19.9678 

218 

0.372 

2/20 

13 

19.9778 

166 

0.375 

2/20 

15 

19.9758 

193 

0.371 

2/21 

15 

19.9837 

192 

0.373 

2/21 

13 

19.9996 

168 

0.370 

2/22 

15 

19.9936 

192 

0.373 

Average  e = 0.370  mcal/vs 



o (Single  measurement)  = 

1 J 

2 X 10  3 = 0.61% 

1 
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Table  A-6.  Long-Term  Calibration  Constants,  0.03  mv  Range  (with  Rotation) 


Date 

t/s 

I/ma 

A/ vs 

e/mcal/vs 

10/18 

13 

9.8995 

135 

0.1 13 

10/22 

17 

9.8856 

178 

0.112 

10/26 

11 

9.9420 

119 

0.110 

10/29 

11 

9.8898 

115 

0.112 

11/19 

17 

9.8266 

176 

0.112 

11/19 

17 

9.8210 

176 

0.112 

11/29 

17 

9.9364 

182 

0.111 

12/3 

20 

9.7795 

207 

0.109 

12/3 

22 

9.7801 

225 

0.111 

12/3 

25 

9.7825 

257 

0.1 11 

12/4 

13 

9.8216 

137 

0.109 

12/5 

14 

9.9319 

153 

0.108 

12/14 

20 

9.9684 

213 

0.112 

2/8 

19 

9.9909 

205 

0.1 11 

2/21 

20 

10.0048 

214 

0.112 

2/28 

16 

10.0068 

171 

0.112 

3/1 

15 

10.0068 

162 

0.111 

Average  e = 0. 1 1 1 mcal/vs 

a (5 

ingle  measurement)  = 1 

X 

o 

1 

u> 

II 

.3% 

Table  A-7.  Long-Term  Calibration  Constants,  0.01  mv  Range  (with  Rotation) 


Date 

t/s 

I/ma 

A/vs 

e/mcal/vs 

10/18 

13 

5.9957 

152 

0.0367 

11/1 

12 

5.9764 

137 

0.0375 

11/26 

14 

5.9959 

165 

0.0366 

1 1/29 

18 

5.9933 

208 

0.0372 

1/30 

14 

5.9949 

156 

0.0387 

3/1 

15 

5.9826 

172 

0.0373 

3/1 

15 

5.9727 

171 

0.0373 

Average  e = 0.0373  mcal/vs 

o ( 

Single  measurement)  = 7 

X 10'4=  : 

.8% 
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Table  A-8.  Keithley  Range  Calibration  Constants 

i 


Range 

eexpta 

Range  factor*5 

Range  c 

ecalc  ^ 

mv 

mcal/vs 

mv 

mcal/vs 

10 

36.971  (30) e 

1.0 

10 

36.971 

3 

11.088  (11) 

1.0032  (l)f 

2.990 

11.056 

1 

3.697  ( 7) 

1.0038  (2) 

0.996 

3.683 

0.3 

1.112  ( 6) 

0.9999  (2) 

0.300 

1.109 

0.1 

0.370  ( 2) 

0.999  (3) 

0.100 

0.370 

0.03 

0.111  ( 1) 

0.997  (1) 

0.301 

0.111 

0.01 

0.0373  ( 7) 

1.0 

0.01 

0.0370 

a Average  calibration  constants  from  long-term  calibration  experiments  shown  in  tables  4-10. 
b Ratio  of  indicated  calorimeter  full-scale  voltage  output  to  actual  ESI-measured  voltage  - average  or 


two  experiments. 

c Adjusted-Keithley  ranges  obtained  by  dividing  the  original  Keithley  ranges  by  range  factors. 

^ Calculated  calibration  constant  obtained  by  using  the  average  10  mv  experimentally  determined  calibration 
constant  in  equation  8. 

e The  numbers  in  the  parentheses  are  the  standard  deviations  of  the  last  digit  for  a single  measurement,  i.e. . 
36.971  (30)  = 36.941  to  37.001. 

r 

The  numbers  in  the  parentheses  are  the  experimental  spread  in  the  last  digit  for  two  measurements,  i.e. . 
1.0032(1)=  1.0031  to  1.0033. 
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Table  A-9.  Area  Averages  for  Calibration  Curves,  1-pv  Range  (with  Rotation) 


Curve  section  time  a 

Average  curve  areas  b 

Class  A c 

Class  B d 

min 

vs 

2-15-5 

144  ( 6)  e 

114(14) 

2-10-4 

138  ( 2) 

105  ( 5) 

2-8-4 

133  ( 1) 

103  ( 6) 

2-6-4 

121  ( 2) 

93  ( 9) 

2-10-2 

137 ( 2) 

106  ( 9) 

2-8-2 

131  ( 1) 

101  ( 6) 

2-6-2 

1 16  ( 1) 

91  ( 5) 

0-10-2 

110(11) 

102  (31) 

0-8-2 

106(17) 

61  (34) 

0-6-2 

64(15) 

46  ( 7) 

a The  areas  listed  for  each  curve  section  time  were  computed  from  the  same  experimental 
curves. 

11  Based  on  the  average  10-mv  calibration  constant,  the  area  that  should  be  observed  if  there 
is  no  heat  of  friction  or  electrical  noise  is  146  volt  seconds. 

c Area  averages  of  six  heat  curves  characterized  by  figure  3.A. 

d Area  averages  of  four  heat  curves  characterized  by  figure  3,B. 

L The  numbers  in  the  parentheses  are  the  standard  deviations  of  a single  area  measurement, 
z.e.,  144(6)=  138  to  150. 


Table  A-10.  Calorimeter  Sensitivity,  1-pv  Range  (with  Rotation) 


Experiment 

No. 

I/ma 

t/s 

Aa/vs 

%bsb 

qobs-qavg 

peal 

peal 

1 

0.9999 

35 

mm 

377 

20 

2 

0.9999 

35 

WL'M 

377 

20 

3 

0.9999 

35 

95 

352 

-5 

4 

0.9999 

35 

94 

348 

-9 

5 

0.9999 

35 

89 

329 

-28 

Average  q = 

357  peal 

Sensitivity  = o (Si 

ngle  measi 

jrement  in 

q)  = ±21 

/-teal 

a 2-15-5  Curve  section  time;e^  = 3.697  peal, 
k Electrical  heat  added  was  419  peal. 
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Tabic  A-l  1.  Experimental  Heats  of  Dilution 
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Experimental  error. 


Table  A- 1 1 . (Contd) 
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Table  A-l  I . (Con til) 
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‘Data  point  outside  2 a in  least  squares  fit. 
‘Wet  wall  experiment, 
t Experimental  error. 
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Experimental  error. 


Table  A- 1 1 . (Contd) 
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outside  2 a in  least  squares  fit. 


Table  A- 1 1 . (Contd) 
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0.5023  1.100  7.8  1.8  0.1709  26.62  26.84  26.72 

2.0604  3.820  0 0 0.6251  104.40  104.40  104.40 
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Wet  wall  experiment. 


T able  A- 1 1 . (Contd) 


Experimental  error. 
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0.4303  1.012  6.7  1.6  0.1524  32.97  32.99  33.09 

0.1524  0.379  2.5  0.6  0.0559  10.67  11.63  10.72 

2.2847  4.343  0 0 0.7082  168.67  168.50  168.67 

0.7082  1.567  10.8  2.5  0.2393  54.60  54.82  54.80 

0.2393  0.585  3.8  0.9  0.0867  17.72  18.31  17.79 
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0.5170  2.009  4.1  1.8  0.3255  59.91  60.10  60.10 

1.0738  4.134  0 0 0.6746  143.60  140.76  143.60 

0.6746  2.548  5.3  2.4  0.4259  84.25  81.47  84.51 

0.4259  0.898  6.7  1.5  0.1386  86.16  85.20  86.43 
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Data  point  outside  2 a in  least  squares  fit. 


(Contd) 
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Wet  wall  experiment. 

Data  point  outside  2 a in  least  squares  fit. 


Table  A- 1 1 . (Contd) 
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Table  A-ll.  (Contd) 
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2.6869  4.427  0 0 0.7368  177.81  177.44  177.81 

0.7368  1.682  10.2  2.3  0.2497  51.06  50.77  51.22 

0.2497  0.676  3.9  0.9  0.0976  16.14  16.39  16.20 


Table  A-14.  Pairwise  Enthalpies  of  Interaction  {ab}^ 


Compound 

|ab}h  calculated 

|ab|  h experimental 

a 

b 

cal  kg/mol2 

cal  kg/mol  2 

Methanol 

Methanol 

102 

591 

Ethanol 

Ethanol 

128 

591,582 

n-Propanol 

n-Propanol 

169 

121 1 , 1342 

n-Butanol 

n-Butanol 

225 

2803, 2402 

t-  Butanol 

t-Butanol 

225 

1 572,  1574,  145 5 

Ethanol 

n-Propanol 

144 

892 

Ethanol 

n-Butanol 

161 

1112 

Ethanol 

t-  Butanol 

161 

882 

n-Propanol 

n-Butanol 

202 

1792 

n-Propanol 

t-  Butanol 

202 

1582 

Isopentanol 

Isopentanol 

297 

4203 

Glycerol 

Glycerol 

104 

603 

Glycolamide 

Glycolamide 

-10 

-466 

1 Reference  1. 


Franks,  F.  Private  communication,  1975. 

3 Reference  2. 

4 Reference  3. 

3 Reference  4. 

^ Reference  5. 
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Table  A-l  S.  Magnitude  of  Terms  in  Equation  1 2 

SBAA/T  _ JAA[hV;  - «RT2{AAjgvj  _ rt2^e(A) 
6 1/T  RT 


Compound 

|AA|g 

<*RT2{AA}gV; 

RT20e(A) 

RT 

cal  8/mol2 

cal/mol2 

cal  8/mol 2 

cal  8/mol2 

Urea 

-83.9 

-26 

-2.1 

18 

Sucrose 

138.1 

42 

3.3 

36 

Urea-sucrose 

-143.4 

-37 

-3.0 

27 

Table  A-l  6.  Pairwise  Enthalpies  of  Interaction 


R-R 

+ 

R-A 

+ 

> 

1 

> 

- 

R-OH 

+ 

OH -OH 

+ 

A-OH 

_ 

R = hydrocarbon  group  - methyl  or  methylene. 
A = amide  group. 

OH  = hydroxy  group. 
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Table  A-17.  Experimental  jab}^  exp  Versus  Calculated  jab  ^ calc 
Pairwise  Heats  of  Interaction 


Solute  a-b 

jabj^  calc 

|ab}h  exp 

cal  kg/mol2 

cal  kg/mol2 

NMF-NMF 

50 

65 

NMF-NMA 

72 

88 

NMF-NMP 

94 

129 

NMF-NBA 

137 

211 

NMA-NMA 

101 

56 

NMA-NMP 

131 

94 

NMA-NBA 

189 

150 

NMP-NMP 

168 

152 

NMP-NBA 

242 

258 

NBA-NBA 

346 

353 

U-U 

-96 

-84 

NMF-U 

-9 

-26 

NMA-U 

13 

0 

NMP-U 

34 

43 

NBA-U 

77 

63 

EG-EG 

103 

86 

EG-P 

104 

103 

EG-G 

105 

102 

EG-S 

108 

164 

P-P 

106 

94 

P-G 

108 

110 

P-S 

115 

150 

G-G 

111 

82 

G-S 

127 

113 

S-S 

140 

138 

EG-U 

64 

11 

P-U 

27 

-26 

G-U 

-10 

-90 

S-U 

-122 

-143 

EG-NMF 

95 

99 

EG-NMA 

114 

125 

EG-NMP 

132 

168 

EG-NBA 

168 

203 

* These  values  were  calculated  by  substituting  the  following  least  square  parameters  into 
equation  (19):  CQ  = 102,  EcH2-CH2  = 7'6’  EA-A  = '88>  ECHOH-CHOH  = 0 3> 


\ ECH2-A  = ,4>  ECH2-CHOH  = g ECHOH-A  = 12 

NMF  = N-methylformamide,  NMA  = N-methylacetamide,  NMP  = N-methylpropionamidc, 
; NBA  = N-butylacetamide,  U = urea. 

1 EG  = Ethylene  glycol,  P = Pentaerythritol,  G = Glucose,  S = Sucrose. 
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Table  A-17.  (Contd) 


Solute  a-b 

jab}^  calc 

jab(h  exp 

cal  kg/mol^ 

cal  kg  mol^ 

P-NMF 

89 

88 

P-NMA 

125 

118 

P-NMP 

162 

137 

P-NBA 

234 

159 

G-NMF 

82 

95 

G-NMA 

137 

176 

G-NMP 

191 

251 

G-NBA 

301 

309 

S-NMF 

62 

12) 

S-NMA 

172 

237 

S-NMP 

281 

327 

S-NBA 

500 

424 
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COMPUTER  PROGRAM  LISTING 


The  program  listing  contained  in  this  appendix  was  used  to  obtain  the  areas  in 
volt-seconds  under  the  various  heat  of  dilution  curves.  The  computer  used  to  run  the  program 
was  a Data  General  Nova  1200.  The  calorimeter  output  for  each  dilution  was  stored  on  punched 
paper  tape  which  was  then  read  into  the  computer  and  stored  as  a data  file  on  a magnetic  disk. 
The  program  is  set  up  to  least-square  the  heat  curve  baseline  and  use  Simpson’s  rule  to  calculate 
the  area  under  the  curve.  These  areas  would  then  be  multiplied  by  the  appropriate  calibration 
constant  to  obtain  the  heat  of  a dilution. 


The  data  file  consists  of  a set  of  voltage  values  listed  in  five  columns,  each  column 
containing  six  digit  numbers  with  a sign  and  decimal  point.  The  last  digit  of  each  number  serves 
to  signal  the  computer  that  a dilution  with  mixing  was  initiated  when  this  digit  is  unity.  At  all 
other  times,  the  last  digit  is  zero.  An  abbreviated  example  of  a data  file  is  shown  below: 


+0.00100 

+0.00120 

+0.00101 


+0.00120 

+0.00110 

+0.00900 


+0.00100 

+0.00120 

+0.01800 


+0.00120 

+0.00110 

+0.02700 


+0.00100 

+0.00100 

+0.00360 

etc. 


The  program  itself  has  been  written  in  such  a way  as  to  require  no  further  elaboration 
since  descriptive  comments  have  been  included  throughout. 


Program  Listing 


LKB  LINEAR  LEAST  SQUARE  AND  SIMPSONS  RULE  PROGRAM 
THIS  PROGRAM  READS  IN  EXPERIMENTAL  VOLTAGES,  DOES 
A LEAST  SQUARES  ON  THE  EASELINE,  CALCULATES  THE 
AREA  UNDER  THE  BASELINE  PORTION  OF  THE  CURVE  AND 
SUBTRACTS  THIS  AREA  FROM  THE  AREA  UNDER  THE  CURVE 
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D I Mfc  NS I ON  V ( 600 ) , KM ( 600 ) , I DATA (11) 

DET  < Tl. T2( T3; T4 > =T 1*T4-T2*T3 

ACCEPT  " THE  NUMBER  OF  POINTS  BEFORE  MIXING  K=  ", K, 

"THE  NUMBER  OF  INTERVALS  UNDER  THE  CURVE  N(EVEN)=  ' • N, 

"THE  NUMBER  OF  POINTS  AFTER  THE  CURVE  BEFORE  LEASTd5> 
SQUARING  BEGINS  NN=  ",NN, "THE  NUMBER  OF  POINTS  AFTERS 1 5 
THE  CURVE  USED  IN  THE  LEAST  SQUARES  IS  M=  ",  M, 

TIME  INTERVAL  BETWEEN  EACH  POINT  IT I ME 
" I TEST * ",  I TEST 
TYPE  "INPUT  DATA  FILE  NAME" 

READ (11, 100)  I DATA ( 1 ) 

(■  . - 1 i T ( S 2 O ) 

A i 1 FOPEN  (1,1 DAT A) 

RFWTND  1 

L ~ i.  . 

N 1 =L+4 

READd  101 , END=20 ) < V ( I ) , KM < I ) , I=L,  N1  ) 

FORMAT  (F7  4,  1 1 , 4 ( F9  4,  1 1 ) > 
v I TEST  EQ.  0)  GO  TO  7 
, i'TEUO, 106)  (V(  I >,  KM(  I >,  I=L,  Nl) 

< IX  5 ( F 1 0 5,  2X,  1 1 ) ) 

L-L+5 
GO  TO  10 

ALL  FCLOS  ( 1 ) 

U r FOR  START  OF  EXPERIMENT 

,_l~0 
J=.J+ 1 

IF  < J GE  L)  GO  TO  200 
iFt'KM(J)  EQ.  0)  GO  TO  22 
WRITE( 10, 102)  J 

FORMAT < //IX, "THE  START  OF  MIXING  OCCURED  AT  POINT  J =",1103 

l IMV=0.  0 
3UMT— 0.  0 
SI  IMTT =0.  O 
SUMVT=0  0 
SI.IMN-O.  0 


r t-ouare  the  baseline 


N2=J— K 

N3= J+M+N+ 1 +NN 
DO  40  I=N2,  N3,  1 
T=I 

N4= J+N+ 1 +NN 

IF ( I.  GE.  J.  AND.  I.  LE.  N4)  GO  TO  40 
SUMV=SUMV  + V(I) 

SUMT  = SIJMT  + T 
SUMTT  = SUMTT  + T*T 
SUMVT  = SUMVT  + T*V(I) 

SUMN=SUMN+1.  0 
40  CONTINUE 

CRAMERS  RULE  PERFORMED  ON 

A ( K+M ) + B(SUMT)  = SUMV 

A (SUMT)  + B( SUMTT)  = SUMVT 

THE  ABOVE  DETERMINANT  IS  DERIVED  FROM 

V( I ) = A + B*T< I ) 

All  = SUMN 
A22  = SUMTT 
A 12  = SUMT 
A21  = SUMT 
Cl  = SUMV 
C2  = SUMVT 

DET1  = DET (All,  A12*  A21<  A22 ) 

DET2  = DET(C1,  A12,  C2,  A22) 

DET3  = DET(AU,  Cl,  A12,  C2) 

A = DET2/DET 1 
B = DET3/DET 1 
WRITE( 10,  103)  A,  B 

103  FORMAT (//IX,  "THE  LEAST  SQUARE  CONSTANTS  ARE:  ", 
1//1X,  "A  =",  F15.  9,  //IX,  "B  =" , F15.  9) 

FN=N 

FJ=J 

AREAB  = A*FN  +B*FJ*FN  + ( B/2 ) * ( FN*FN ) 

TIME=ITIME 

AREAB=AREAB*T I ME 

IF (I  TEST  EQ.  1 ) TYPE " AREAB " , AREAB 
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DO  SIMPSONS  RULE  ON  CURVE  VOLTAGE 

SUM  = V(J) 

N5=J+1 
N6= J+N-3 

DO  50  I = N5,  N6,  2 
50  SUM  = SUM+  4*V( I ) +2*V( 1+1 ) 

SUM  = SUM+4*V( J+N-l )+V( J+N) 

IF<  I TEST.  EQ.  1 ) TYPE  "SUM" , SUM 

ARE AC  = (TIME/3.  0)*SUM 

IF  ( I TEST.  EQ.  1 ) TYPE  "AREAC",  AREAC 

AREA  = AREAC-AREAB 

IFdTEST.  EQ.  1)  TYPE  "AREA",  AREA 

WRITE( 10, 104)  AREA.  AREAC,  AREAB 

104  FORMAT (//IX,  "AREA  = ",  F15.  9, //IX,  "AREAC  = ",F15.  9, 
1//1X,  "AREAB  = " , F15.  9) 

IF(  ITEST.  EQ.  0)  GO  TO  70 

WR I TE ( 10,  105)  K,  N,  M,  All,  A12,  A21,  A22,  Cl,  C2 

105  FORMAT  < //IX , "K=",  14,  /IX,  "N=",  15,  /IX,  "M=",  15, 

2 /IX,  "All  = ",  F 10.  2.  /IX,  "A12=" , F10.  1,  /IX,  "A21  = ",  F10.  1, 

3 /IX,  "A22=",  FI 5.  2,  /IX,  "Cl  = "»  F15.  7,  /IX,  "C2=",  F20.  7) 
70  CONTINUE 

SEARCH  FOR  NEXT  START  OF  EXPERIMENT 

J=.J+10 

GO  TO  22 

200  ACCEPT  "NEW  TIME  INTERVAL  NT,  IF  ZERO  STOP,  = ",  NT 
IF  (NT.  EQ.  0)  GO  TO  300 
DO  210  I =NT , 3000,  NT 
IN=I/NT 
V( IN ) =V ( I ) 

210  KM (IN) =KM ( I ) 

I T I ME=NT  * I T I ME 
K=K/NT 
M=M/NT 
NN=NN/NT 
N=N/NT 
L=L/NT 
GO  TO  20 
300  STOF 
END 
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DISTRIBUTION  LIST  FOR  ARCSL-SP-77003 


Names 

Copies 

Names 

CHEMICAL  SYSTEMS  LABORATORY 

OFFICE  OF  THE  DIRECTOR 
Attn:  DRDAR-CLP 

1 

US  ARMY  TEST  & EVALUATION  COMMAND 
Record  Copy 

CB  DETECTION  & ALARMS  DIVISION 

CDR,  APG 

Attn:  STEAP-AD-R/RHA 

Attn:  DRDAR-CLC 

1 

APG-Edgewood  Area,  Bldg  E5179 

DEVELOPMENTAL  SUPPORT  DIVISION 

Commander 

Attn:  DRDAR-CLJ-I 

4 

US  Army  Test  & Evaluation  Command 

ENVIRONMENTAL  TECHNOLOGY  DIVISION 

Attn.  STEAP-CP-T 
APG-Aberdeen  Area 

Attn:  DRDAR-CLT 

1 

MUNITIONS  DIVISION 

Commander 

Dugway  Proving  Ground 

Attn:  DRDAR-CLN 

1 

1 

Attn.  STEDP-PO 

PHYSICAL  PROTECTION  DIVISION 
Attn:  DRDAR-CLW 

Dugway,  UT  84022 
DEPARTMENT  OF  THE  NAVY 

RESEARCH  DIVISION 

Chief  of  Naval  Research 

Attn:  DRDAR-CLB 

I 

Attn:  Code  443 

Attn:  DRDAR-CLB-B 

I 

800  N.  Quincy  Street 

Attn:  DRDAR-CLB-C 

1 

Arlington,  VA  22217 

Attn:  DRDAR-CLB-P 

1 

Attn:  DRDAR-CLB-P  (Author's  copies) 

10 

Commander 

Attn:  DRDAR-CLB-R 

I 

Naval  Facilities  Engineering  Command 

Attn:  DRDAR-CLB-T 

1 

Attn:  Code  03 

DEPARTMENT  OF  DEFENSE 

200  Stovall  Street 
Alexandria,  VA  22332 

Administrator 

Defense  Documentation  Center 

DEPARTMENT  OF  THE  AIR  FORCE 

Attn:  Accessions  Division 

12 

Cameron  Station 

HQ,  USAF/SGPR 

Alexandria,  VA  22314 

Forrestal  Bldg 

US  ARMY  MATERIEL  DEVELOPMENT  AND 
READINESS  COMMAND 

Commander 

US  Army  Materiel  Development  and 
Readiness  Command 
Attn:  DRCSF-S 

1 

WASH  DC  20314 

5001  Eisenhower  Ave 
Alexandria,  VA  22333 

US  ARMY  ARMAMENT  RESEARCH  AND 
DEVELOPMENT  COMMAND 

Commander 

US  Army  Aberdeen  Proving  Ground 
Attn:  DRDARQA 

1 

Aberdeen  Proving  Ground,  MD  21010 
Director 

Ballistic  Research  Laboratory 

Attn:  DRDAR-TB  (Dr.  David  Hansen) 

1 

Building  328 

Aberdeen  Proving  Ground,  MD  21005 


Copies 


1 


1 


2 


I 


1 


1 


Jd 


89 


